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ABSTRACT 


Two  dimensional  steady  symmetric  incompressible  laminar  flow  past 
a  class  of  blunt  and  sharp  nosed  bodies  is  investigated  in  optimal  co¬ 
ordinates.  The  analysis  is  carried  out  for  different  problem  parameters 
and  the  solution  is  specialized  for  the  cases  of  the  parabola,  a  semi- 
infinite  thin  flat  plate  and  the  flow  against  a  vertical  wall. 

The  problem  is  formulated  by  mapping  the  body  from  a  Cartesian 
plane  into  a  conformal  plane  by  applying  a  Schwarz-Christoffel  trans¬ 
formation.  <This  type  of  transformation  maps  the  flow  field  into  a 
stagnation  point  flow  in  the  conformal  plane.  The  governing  differential 
equations  are  developed  in  an  optimal  coordinate  system  such  that  the 
viscous  solutions  contain  the  first-order  inviscid  solution  as  well  as 
the  solution  for  l\bv  due  to  displacement  thickness.  This  type  of 
formulation  of  the  viscous  flow  problem  in  optimal  coordinates  results 
in  a  coupled^Set  of  flow  and  coordinate  equations.  A  boundary-layer 
and  pardholized  Navier-Stokes  analysis  is  used  for  the  solution  procedure 
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i^the  present  optimal  coordinate  system. 

Optimal  coordinates  are  computed  according  to  the  classical  first- 
order  boundary-layer  approximation  as  well  as  with  a  parabolized  version 


of  streamfunction  vorticity  form  of  the  full  Navier-Stokes  equations. 

A  boundary-layer  type  forward  marching  numerical  scheme  is  employ^d^to 


solve  the  flow  equations  in  optimal  coordinates. 
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A  The  analysis  is  carried  out  for  two  example  problems,  a  semi-infinite 
thick  plate  and  a  semi-infinite  blunted  wedge.  The  solutions  are  obtained 
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.  for  different  body  geometries  (bluntness  parameters)  in  both  examples. 
Results  for  skin  friction,  displacement  thickness,  pressure  gradient 
parameter  and  optimal  coordinates,  for  different  problem  parameters, 
are  presented  for  unseparated  flow  cases.  Physical  quantities  such  as 
surface  pressure  gradient  distribution,  skin  friction  and  displacement 


thickness  for  cases  of  flow  past  the  parabola,  the  thin  flat  plate  and 
flow  against  a  vertical  wall  are  compared  with  existing  numerical  and 
analytical  results.  _ 
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I.  INTRODUCTION 


A  suitable  choice  of  a  coordinate  system  can  simplify  the  complexity 
of  a  problem  especially  one  involving  non-linear  partial  differential 
equations.  In  some  cases  the  proper  choice  of  independent  and  dependent 
variables  reduces  the  original  partial  differential  equation  to  an  ordi¬ 
nary  differential  equation,  i.e.  if  similarity  conditions  exist.  The 
Falkner-Skan  similarity  variables  are  an  example  of  such  a  case  where 
the  boundary- layer  equations  reduce  to  an  ordinary  differential  equation 
called  the  Falkner-Skan  equation.  Due  to  the  simplicity  of  ordinary 
differential  equations,  the  numerical  solution  of  ordinary  differential 
equations  can  be  found  more  accurately  and  with  less  effort,  in  general, 
than  the  numerical  solution  of  the  original  partial  differential  equations. 

In  many  singular  perturbation  problems,  which  show  a  boundary-layer 
like  behavior,  the  coordinate  system  chosen  for  the  problem  may  play  an 
important  role  in  describing  the  flow  appropriately.  In  other  words,  the 
solution  of  boundary-layer  type  flows  may  depend  on  the  coordinate  system 
chosen.  This  fact  did  not  receive  a  systematic  investigation  until  Kaplun^ 
introduced  the  idea  of  optimal  coordinates  which  according  to  Kaplun 
removes  the  singularity .due  to  the  boundary-layer  approximation  and  renders 
the  boundary-layer  solution  uniformly  valid  in  the  entire  flow  field  for 
high  Reynolds  number  flow. 

The  idea  of  optimal  coordinates,  although  first  mentioned  by  Kaplun  is 

2  3  4  5 

also  mentioned  in  the  works  of  Latta  ,  Friedrichs  ,  Weyl  and  Segel  . 

The  dependence  of  the  boundary-layer  solution  on  the  coordinate  system  is 
discussed  by  Van  Dyfck^  for  flow  past  a  flat  plate  in  Cartesian  and  para¬ 
bolic  coordinates.  Some  other  applications  of  optimal  coordinate  theory 
liave  been  given  by  Crespo  da  Silva  and  Davis^ . 


Kaplun*  shows  that  the  inviscid  high  Reynolds  number  solution  is 
independent  of  the  coordinate  system,  since  it  is  governed  by  a  vector 
equation.  However,  he  also  shows  that  in  a  boundary-layer  type  analysis, 
while  going  from  the  viscous  flow  region  to  the  inviscid  flow,  the 
boundary- layer  solution  produces  different  types  of  outer  inviscid  flows 
depending  on  which  coordinate  system  has  been  chosen.  The  coordinate 
system  which  provides  the  outer  flow  including  the  displacement  effect, 
to  a  given  order,  is  an  optimal  coordinate  system  to  that  order.  There¬ 
fore  Kaplun's  investigation  of  the  dependence  of  the  boundary-layer  solu¬ 
tion  on  the  coordinate  system  leads  to  the  possibility  of  developing  the 
most  appropriate  coordinate  system  for  carrying  out  high  Reynolds  number 
flow  computations. 

Since  the  introduction  of  the  idea  of  optimal  coordinates  by  Kaplun, 
several  investigators  have  explored  the  application  of  optimal  coordinate 
theory  to  obtain  uniformly  valid  solutions  for  singular  perturbation 
problems  in  fluid  mechanics.  A  review  of  the  literature  related  to  optimal 
coordinates  is  presented  by  Legner^.  Legner^  shows  that  the  application 
of  Kaplun's  idea  is  not  limited  to  fluid  mechanics  problems,  but  also  can 
successfully  be  applied  to  singular  perturbation  problems  in  other  applied 

g 

engineering  areas.  For  example,  Zauderer's  work  on  defraction  problems 
involving  multi-region  transition  and  Segel's  work  on  the  vibration  of 
convex  cylinders  in  a  uniform  stream  are  worth  mentioning. 

9 

An  extension  of  Kaplun's  correlation  theorem  was  made  by  Van  Dyke 
who  extended  Kaplun's  first-order  theorem  to  second-order.  This  allows 
one  to  transform  a  second-order  boundary- layer  solution  from  one  coordinate 


2 


system  to  another.  A  further  development  of  Kaplun's  work  was  presented 
by  Legner^  who  extended  the  entire  optimal  coordinate  theory  to  higher 
order.  Subsequent  generalization  of  Legner's  approach  to  finding  optimal 
coordinates  was  provided  by  Davis. ^ 

The  numerical  determination  fo  coordinate  systems  for  use  with  the 

boundary-layer  or  Navier-Stokes  equations  has  been  studied  by  several 
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investigators,  see  for  example,  Thompson  et  al.  ’  ’  In  these  methods 

in  order  to  generate  a  coordinate  system,  a  criterion  must  be  specified 
which  the  coordinate  system  should  satisfy.  Therefore,  if  we  set  our 
criterion  for  the  coordinate  system  such  that  the  boundary-layer  solution 
in  these  sets  of  coordinates  reproduce  the  inviscid  solution  including 
the  displacement  effect,  the  coordinates  thus  obtained  will  be  optimal  in 
the  sense  of  Kaplun's  definition.  Davis^°  has  discussed  how  the  optimal 
coordinates  can  be  generated  numerically  in  a  manner  similar  to  that  used 
by  Thompson  et  al.  ^  for  non  optimal  coordinates. 

The  advantages  of  the  use  of  the  optimal  coordinates  with  the 
boundary-layer  and  the  Navier-Stokes  equations  for  obtaining  numerical 
as  well  as  analytical  solutions  are  ellaborated  in  the  works  of  Davis^ 
and  Davis. ^  In  the  optimal  coordinate  study  by  Davis^  it  has  been 
shown  that  for  a  plane  symmetric  flow,  the  use  of  the  Legner  condition 
and  the  assumption  of  the  orthogonality  condition  for  the  coordinate 
system  results  in  a  set  of  conformal  coordinates.  These  conformal  co¬ 
ordinates,  which  in  the  sense  of  optimal  coordinate  theory  are  the  best 
coordinates,  are  used  in  the  numerical  scheme  of  Davis^  for  viscous 
flow  past  the  parabola,  which  becomes  a  flat  plate  as  the  nose  radius  of 
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curvature  approaches  zero.  This  coordinate  system  shows  considerable 
improvement  over  other  coordinate  systems  in  the  convergence  and  accur¬ 
acy  of  the  numerical  solution.  These  coordinates  were  also  successful 
in  removing  the  leading  edge  singularity  for  the  flow  past  the  flat  plate. 

However,  the  use  of  conformal  coordinates  in  the  case  of  flow  past  wedges 
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(see  Davis,  U.  Ghia  and  K.  Ghia  )  ,  although  showing  good  resolution  and 
convergence  in  most  of  the  flow  field,  fails  to  remove  the  leading  edge 
singularity  except  for  the  special  cases  of  semi-wedge  angle  equal  to 
zero,  (the  flat  plate  case),  and  90°  (the  vertical  wall  case).  Davis^ 
has  further  shown  that  it  may  be  possible  to  generate  a  non-orthogonal 
optimal  coordinate  system  which  removes  the  leading  edge  singularity  for 
the  sharp  wedge  cases,  by  introducing  the  required  arbitrariness  into 
the  generation  of  optimal  coordinates. 

The  present  numerical  work  will  show  that  in  optimal  coordinates, 
parabolization  of  the  Navier-S tokes  equations  renders  a  considerable 
simplification  in  the  numerical  treatment  of. the  Navier-S tokes  equations 
and  at  the  same  time,  that  it  is  possible  to  achieve  good  agreement 
between  solutions  to  the  parabolized  equations  and  the  full  Navier-S tokes 
equations.  Previous  work  on  developing  parabolized  models  for  the  same 
type  of  problems,  but  not  in  a  fully  optimal  coordinate  system,  has  been 
done  by  Davis, ^  D.  Ghia  and  Davis, ^  and  Werle  and  Bernstein.^ 

In  the  present  analysis,  a  parabolic  form  of  the  Navier-Stokes 
equations  is  therefore  developed  in  a  generalized  non-orthogonal  optimal 
coordinate  system.  The  purpose  of  this  investigation  is  to  determine 
if  there  is  an  advantage  in  using  these  optimal  coordinates. 


Analytical  determination  of  the  present  optimal  coordinate  system 
is  rather  difficulty  since  in  order  to  determine  the  present  set  of 
coordinates  it  is  essential  to  solve  for  the  flow  due  to  displacement 
thickness.  Only  for  a  few  simple  cases  of  viscous  flow  is  the  analytical 
expression  for  the  inviscid  flow  due  to  displacement  thickness  known. 

In  most  cases  the  solution  to  the  inviscid  flow  due  to  displacement  thick¬ 
ness,  must  be  generated  numerically.  Therefore  in  the  present  analysis, 
a  method  is  developed  to  numerically  generate  the  optimal  coordinates 
using  a  finite  difference  technique.  The  flow  equations  are  solved 
utilizing  a  boundary-layer  like  forward  marching  technique  and  the 
coordinate  equation  is  solved  using  a  time  dependent  relaxation  scheme. 

Two  model  problems  are  considered  in  order  to  demonstrate  the  application 
of  the  present  analysis.  It  is  shown  through  model  problems  that  it  is 
possible  to  numerically  generate  a  set  of  optimal  coordinates  and  that 
the  solution  obtained  from  the  present  calculations  shows  good  agreement 
with  the  existing  analytical  and  numerical  results. 
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II.  GENERATION  OF  OPTIMAL  COORDINATES 

It  was  mentioned  in  the  Introduction  chat  in  order  to  determine 
optimal  coordinates  in  a  high  Reynolds  number  flow  problem,  it  is 
necessary  to  know  apriori  not  only  the  first-order  inviscid  flow  but 
also  the  flow  due  to  displacement  thickness.  This  implies  that  one 
should  have  the  first-order  boundary-layer  solution  at  hand  in  order  to 
determine  optimal  coordinates  to  first  order.  This  information  is 
available  in  the  form  of  self-similar  solutions  for  a  few  simple  flows 
but  in  most  cases  one  must  resort  to  a  numerical  method  in  order  to  de¬ 
termine  the  boundary-layer  solution.  The  inviscid  flow  due  to  the  dis¬ 
placement  effect  is  then  determined  from  a  separate  analysis  or  numerical 
calculation. 

In  the  subsequent  section  a  brief  description  is  presented  on  how 
to  generate  optimal  coordinates  numerically  for  a  two  dimensional  in¬ 
compressible  symmetric  flow.  However,  the  analysis  could  be  extended 
to  unsymmetric,  compressible,  three  dimensional,  etc.  flows  as  well, 
see  Legner.^ 

For  the  present  study  we  will  assume  that  the  form  of  the  optimal 
coordinates  chosen  are  such  that  they  produce  a  separable  solution  for 
the  inviscid  flow  to  arbitrary  order.  These  types  of  optimal  coordinates 
are  described  and  developed  in  detail  by  Davis. ^  However,  a  brief 
summary  of  Davis  ^  analysis  is  presented  in  the  following  section  for 
the  sake  of  completeness. 

In  the  Davis^  analysis  it  is  assumed  that  the  Independent  variables 
are  chosen  such  that  a  separable  solution  exists  in  the  inviscid  while 
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an  unseparable  solution  is  confined  at  most  to  the  viscous  region.  In 
some  special  cases  the  viscous  region  may  also  be  separable.  Choice  of 
such  variables  represent  definite  advantages  in  determining  solutions 
analytically  as  well  as  numerically. 

In  order  to  obtain  a  set  of  optimal  coordinates,  hereafter  denoted 
by  a  (5,n)  coordinate  system,  the  flow  field  in  the  physical  (x,y)  plane 
is  mapped  into  a  stagnation  point  type  flow  by  a  conformal  transformation 
of  the  following  form  (see  Figure  1) : 

2  *  f (p+iq)  «  f(^)  .  (2.1) 

The  quantity  5  *  p+iq  represents  the  conformal  plane  and  f  is  an  analytic 
function.  This  type  of  mapping  is  used  in  order  to  simplify  the  boundary 
conditions  which  result  from  relating  the  optimal  to  the  conformal  system 
of  coordinates.  In  addition  there  are  certain  other  advantages  to  this 
type  of  mapping  and  these  will  be  elaborated  on  in  a  later  section. 

The  outer  inviscid  flow,  which  according  to  optimal  coordinate 

theory,  (see,  Legner^ )  is  independent  of  the  coordinate  system,  is  matched 

with  the  boundary-layer  solution  in  the  (C,n)  coordinate  system,  by  apply- 

19 

Ing  the  asymptotic -matching  principle,  (see.  Van  Dyke  ) .  In  the  present 
matching  of  two  solutions  it  is  assumed  that  the  boundary-layer  solution 
in  the  (£,n)  optimal  coordinate  system  is  uniformly  valid  far  from  the 
body,  i.e.  in  the  transverse  direction.  This  type  of  matching  of  the 
two  solutions  results  in  a  set  of  conditions  for  determination  of  the 
optimal  coordinate  system,  the  order  of  which  depends  on  the  number  of 
terms  considered  in  the  expansion  in  terms  of  1/i^Re  as  the  perturbation 
parameter.  In  order  to  further  simplify  the  condition  for  the  optimal 
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coordinates  and  to  obtain  a  more  useful  form,  Davis  uses  the  Legner 

form  of  the  outer  inviscid  streamf unction  expansion  and  we  will  adopt 

the  same  procedure  here.  Finally  invoking  the  fact  that  the  inviscid 

solution  is  separable  to  arbitrary  order  in  the  present  set  of  (£,n) 

optimal  coordinates,  the  proper  form  of  the  outer  flow  including  the 

displacement  effect  is  obtained.  Using  this  form  of  the  inviscid  stream- 

function  in  Laplace’s  equation  governing  the  streamf unction,  i.e., 

2 

7  *  0,  together  with  a  suitable  coordinate  transformation  relation 

to  relate  the  conformal  coordinates  p  and  q  with  the  optimal  coordinate 

system  £  and  n,  the  coordinate  equation  is  developed.  Finally  the 

governing  differential  equations  in  terms  of  streamf unction  and  vorticity 

are  developed  in  optimal  coordinates  from  the  corresponding  equations  in 

20 

a  generalized  coordinate  system  as  given  by  Lagerstrom  p.  66,  by 
applying  the  suitable  coordinate  transformation  manipulations.  The 
necessary  coordinate  transformation  relation,  the  coordinate  equations 
to  determine  the  optimal  coordinate  system  and  the  flow  equations  to¬ 
gether  with  the  boundary  conditions  to  solve  these  equations  in  optimal 
coordinates  will  be  presented  in  their  final  forms  as  given  by  Davis. ^ 

The  details  of  the  derivation  will  not  be  repeated  here. 

In  the  present  analysis  the  flow  field  is  mapped  into  a  conformal 
plane  and  the  corresponding  flow  is  then  investigated  in  the  optimal 
plane.  A  typical  example  of  such  a  mapping  of  the  flow  field  is  presented 
in  Figure  1,  in  which  the  two  dimensional  flow  past  a  circular  cylinder 
is  shown  in  the  physical,  aonformal  and  optimal  plane  in  Figure  la,  lb 
and  lc,  respectively  (see  also  Davis^) . 
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The  transformation  of  independent  variables  from  the  conformal  (p,q) 
to  the  optimal  (5,n)  coordinate  system  is  carried  out  using  the  following 
relations 


3_  _  1_  ,  3_  3_v 

3p  “  Jr  %  35  “  \  3rr 


and 


(2.2a,b) 


3  1  ,  3  3  . 

3q  "  Jr  Pn  35  +  P5  3n5  * 

where  is  the  Jacobian  of  the  transformation  matrix  and  is  given  by 


J  ■  p_q  -p  q_ 

r  *5  n  *nM5 


(2.2c) 


The  rest  of  the  coefficients  of  the  metric  tensor  relating  the 
optimal  coordinate  system  (5,n)  to  the  conformal  coordinates  (p,q)  are 
given  as 


°r  *  gr 


2  2 
p„  + 


22 


6r  "  «r12  "  p5pn  +  q5qn 


2  2 
Yr  "  P5  +  q5 


(2.3a,d) 


and 


2  2  2 
J  “  g  ■  a  y  -S  ■  (p„q-p  q_) 
r  ®r  r'r  r  *5  rr  n  5 
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Similar  expressions  can  be  obtained  to  relate  a  Cartesian  coordinate  system 
to  the  optimal  coordinates.  They  are  rewritten  here  as  given  by  Davis^ 

2  2  .2 

a  “  g22  =  xn  +  yn  ’  h  ar  * 


e  -  g 


12 


sn 


X,X  + 

Z  n  5  n 


h2e 


2  2 

V’s 


h2Y. 


(2.4a,e) 


and 


where 


g  -  ay  -  g 


✓  \ 2  *4  2 

(x_y  -x  y_)  >  h  J 

5  n  n  C  r 


h-£i  • 

We  will  now  write  down  the  governing  differential  equations  in  the 
(£,n)  coordinate  system. 

1.  The  Navier-Stokes  Equations 

The  non-dimens ionalization  used  in  the  flow  equations  is  as  follows . 

All  lengths  are  non-dimensionalized  by  the  viscous  length  v/U^,  where  v 

is  the  kinematic  viscosity  and  Ur  is  the  reference  velocity  in  the  flow. 

2 

The  streamf unction  is  non-dimensionalized  by  v  and  the  vorticity  by  U^/v. 

2 

Finally  the  time  is  non-dimensionalized  by  v/U^. 

With  this  non-dimensionalization  the  streamf unction  equation  is  given 
in  the  (£,n)  coordinate  system  as 


26  a  .a 

f - f  +  —  f  +  —  ~  f  »  g 

nn  yr  Cn  Yr  CC  5  Yr  C 


(2.5) 


where  the  streamf unction  ^  is  defined  in  the  spirit  of  similarity  vari¬ 
ables  as  j 

-  5f(€,n)  (2.6) 

and  the  vorticity  iu  is  likewise  defined  as  | 

5Yr  1 

<*>  -  -  -ys  g(5,n)  =  -  5Hg(C,n)  (2.7)  j 

h  J 

r 


where  the  inverse  scale  factor  H  is  given  as 


The  transformed  vorticity  transport  equation  is  given  by 
2H  2g  H  J 


V  H2  ‘'r  H2  «  H  5  Yr  H 


,6  ,86  8  2  J  H 

,  4  /  r \  4  r  ,  r.  ,12  .  r.  r  n 

5  (~}n  “IT"  (~}C  +  72  ^  +  ~  r  <f+£V 

r  r  r  5  r  r 


(2.12) 


J  ^Hr  2a  Hr  2K  Hn 

- £  (1  +  _i)f  +  r — I  _£ - 1  _H - £  f  ]„ 

Yr  u  H  'V8  lYr  H  Yr  H  Yr  V85 


26r  ar  £  “r  1 

'  g€n  +  77  +  I  =  H  gt  * 

The  boundary  conditions  are  given  as  follows.  The  no-slip  and  zero 
injection  conditions  give  at  n  *  0 
f (5,0)  *  0 


(2.13a,b) 


fn(c,o)  -  0  . 


The  asymptotic  boundary  condition  far  from  the  body  gives 
f(5,n)  ->  n-C 


as  n  -*■ 


(2.14a,b) 


g(5,n)  -*■  0 

where  C  is  a  constant  and  from  the  general  definition  of  the  displacement 
thickness  Davis^  has  shown  that  C  is  a  displacement  constant  in  optimal 


coordinates. 
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therefore  is  a  known  quantity. 

It  is  interesting  to  note  that  the  Equations  (2.15)  and  (2.16)  are 

the  similarity  form  of  the  boundary-layer  equations.  This  shows  that 

the  variables  chosen  for  the  present  analysis  are  self-similar  in  nature. 

Also  substituting  f  ■  g  into  the  vorticity  equation  (2.16)  and  inte- 

nn 

grating  the  resulting  equation  with  respect  to  n  once  and  making  use  of 
the  boundary  condition  on  f^  and  the  higher  order  derivatives  of  f  as 
n  •,  the  vorticity  Equation  (2.16)  reduces  to 


13 


(2.18) 


+  ^nn  +  ® 
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which  Is  the  familiar  Falkner-Skan  equation  (see  Schlichting  ). 

Additional  relations  required  to  solve  the  flow  equations  are  given 
in  the  following.  From  Equation  (2.8),  it  can  be  shown  that 


h  Hr 


(2.19) 


_ji  +  _a 

h  Hr 


(2.20) 


This  completes  the  description  of  the  Navier-Stokes  equations 
together  with  the  necessary  boundary  conditions  in  the  (C,n)  coordinate 
system.  The  derivation  of  the  coordinate  equation  will  be  briefly  out¬ 
lined  in  the  subsequent  section. 

2.  The  Coordinate  Equation 

The  coordinate  equation  is  derived  based  on  the  inviscid  flow  in¬ 
formation  in  the  following  way.  It  is  noted  from  the  viscous  stream- 
function  f  given  by  Equation  (2.6)  and  the  outer  edge  condition  on  f 
given  by  Equation  (2.14a)  that  the  outer  inviscid  streamf unction  in 
C,n  coordinates  is  represented  as 


'{’in  ’  C(n-C) 


(2.21) 


If  the  flow  is  irrotational  to  arbitrary  order,  then  the  streamf unction 

2 

should  satisfy  the  Laplace  equation  7  ■  0,  which  can  be  written  in 

our  case  as 


v  (C (n-C))  -  0 


(2.22) 
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Assuming  pq  ■  £n  and  treating  p  in  a  similar  manner  as  we  have  treated 
the  transformed  streamf unction  and  vorticity  w  it  is  possible  to  obtain 

coordinate  relations  of  the  form  given  in  the  following  (see  Davis^) : 
p  -  5t(5.n)  (2.23) 

and 

n  •  (2.24) 

q  *  t(5,n) 

Making  use  of  these  coordinate  relations  (2.23)  and  (2.24)  and  substitut- 

2 

ing  into  (2.22)  and  writing  the  Laplacian  V  in  the  £  and  n  coordinate 
system  (with  the  help  of  the  transformation  of  independent  variable  re¬ 
lations  (2.2a,b)  from  the  conformal  (p,q)  to  the  optimal  (£,ri)  coordinate 
system)  it  is  now  possible  to  show  that  the  coordinate  Equation  (2.22) 
reduces  to 


26  a  _  a 

t  t.  +-£•  t„  +|-E  t  -  o 

nn  Yr  5n  Yr  55  5  Yr  5 


(2.25) 


where  t(£,n)  is  the  coordinate  relation  function  relating  the  two  co¬ 
ordinate  systems  through  the  expressions  given  in  (2.23)  and  (2.24). 

Determination  of  the  (£,n)  coordinate  system  amounts  to  finding  the 
solution  of  the  partial  differential  Equation  (2.25)  with  suitable  boundary 
conditions. 

The  coefficients  af>  6r»  Yr  and  the  Jacobian  J r>  appearing  in  the 
flow  Equations  (2.5,  2.12)  and  the  coordinate  Equation  (2.25)  can  now 
easily  be  expressed  as  functions  of  £  and  n  by  substituting  the  coordinate 
relations  (2.23)  and  (2.24)  into  Equations  (2.3a-d).  This  results  in 

„  .  ,, 2  +_2U2  2ntn  ,  1 

t  t  t 
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(52  +^)t5tn  tc  +  ettn 


2 

+  \)t2  +  2ctt  +  t2 

t  4  4 


(2,26a-d) 


and 


The  boundary  conditions  for  the  coordinate  Equation  (2.25)  are  pres¬ 
cribed  as  follows.  The  coordinate  Equation  (2.25)  is  elliptic  in  nature 
and  therefore  needs  boundary  conditions  around  the  boundaries  in  order  to 
make  the  coordinate  equation  problem  well  posed.  Since  we  will  be  con¬ 
sidering  symmetric  bodies,  the  coordinate  system  is  symmetric  and  therefore 
t^(o.n)  -  o  (2.27) 

The  asymptotic  boundary  condition  in  the  streamwise  and  transverse  direction 
are  given  as 

t(£,n)  1  as  £  -*■  « 

and  (2.28a,b) 

t(5,n)  ^1  as  n  -*■  ” 

The  two  boundary  conditions  given  by  Equations  (2.28a,b)  imply  that  far 
from  the  body  in  either  direction  the  optimal  coordinates  become  conformal. 

The  remaining  one  condition  is  the  boundary  condition  on  t(£,n)  at 
the  surface.  It  was  shown  by  Davis^  that  this  condition  depends  on  the 
displacement  surface  over  the  body  and  hence  it  is  treated  as  an  unknown. 

In  the  following  section  it  will  be  shown  how  an  interacting  type  analysis  is 
necessary  in  order  to  determine  this  missing  boundary  condition  (t(£,0)) 

and  therefore  to  solve  the  flow  equations  in  optimal  coordinates. 
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Referring  back  to  the  boundary  condition  (2.14a)  on  the  streamfunction 


i.e.  f  +  ri  -  C  as  n  +  it  is  generally  treated  as  f  -*•  1  as  n  -*•  00  and 
C  is  determined  as  C(£)  as  n  -*■  »  if  £  and  n  are  not  optimal  coordinates. 
Therefore  under  normal  conditions,  one  has  no  control  over  the  function 
C(£).  However,  it  was  already  mentioned  that  there  was  one  condition 
missing  from  the  coordinate  equation  in  order  to  describe  the  system  in 
closed  form.  In  fact  the  additional  condition  C(£)  -*•  C  as  n  -*■  necessary 
for  determining  optimal  coordinates, is  used  in  determining  the  surface 
boundary  condition  t(£,0).  This  clearly  indicates  that  in  order  to  solve 
the  coordinate  equation  in  closed  form  it  is  essential  to  solve  at  least 
the  first-order  boundary-layer  flow  such  that  the  displacement  surface 
function  C(£)  +  C  as  ri  +  ®.  Since  the  coefficients  a^,  and  etc. 

are  functions  of  t(£,n),  they  are  therefore  unknown  and  are  determined  as 
a  part  of  the  solution.  This  Implies  that  the  flow  equations  (2.5  and  2.12) 
through  these  coefficients  are  coupled  to  the  coordinate  equations  and  the 
coordinate  Equation  (2.25)  is  coupled  to  the  flow  equations  through  •• 
the  boundary  condition  t(£,0),  the  determination  of  which  depends  on  the 
successful  satisfaction  of  the  condition  C(£)  ■+  C  on  the  flow  equation 
solution  as  n  -*•  00 .  Therefore  the  condition  on  the  displacement  surface 
C(£)  -*■  C  serves  as  an  additional  condition  on  the  flow  and  coordinate 
equations  such  that  the  coordinates  are  optimal. 

It  is  now  clear  that  the  streamfunction  equation  boundary  conditions 
are  not  over  specified  if  we  wish  to  prescribe  f  -*  1  and  f  -*•  rj  -  C  as 
n  *►  “  to  determine  optimal  coordinates.  The  extra  boundary  condition  on 
the  streamfunction  f(£,q)  is  replaced  by  the  missing  boundary  condition 
on  t(£,n) .  Therefore  the  Equations  (2.1)  through  (2.28)  represent  a 
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complete  set  of  governing  differential  equations  together  with  the 
necessary  boundary  conditions  in  order  to  determine  the  optimal  coordinates 
for  two-dimensional  plane  symmetric  flow. 

It  is  interesting  to  note  that  the  coordinate  equation  is  formulated 
based  on  the  inviscid  flow  information  whereas  the  viscous  effects  are 
confined  to  the  flow  equations  only.  Because  of  the  inherent  coupling 
between  the  coordinate  equation  in  optimal  coordinates  and  the  flow  equa¬ 
tions,  an  interacting  type  phenomenon  takes  place  while  determining  the 
viscous  flow  solution  in  optimal  coordinates. 

This  brief  description  of  the  development  of  the  optimal  coordinates 
in  the  case  of  two  dimensional  laminar  incompressible  symmetric  flow 
follows  the  more  detailed  analysis  given  by  Davis'^. 

The  development  of  the  parabolized  vorticity  and  the  streamf unction 
equations  and  also  the  boundary-layer  equations  are  presented  in  the 
following . 

3 .  The  Parabolized  Vorticity  and  Stream  Function  Equations 

The  parabolized  vorticity  and  stream  funcation  equations  are  obtained 

from  the  full  Navier-Stokes  equation  by  neglecting  all  the  terms  in  the 

governing  differential  Equations  (2.5)  and  (2.12),  which  include  second 

derivatives  in  the  streamwise  (C)  direction.  Also  the  cross  derivative 

a 

terms  f^  and  g^  and  —  g^  are  neglected.  The  relevant  advantages  of 
the  use  of  the  parabolized  form  of  the  Navier-Stokes  equations  are 
separately  ellaborated  in  the  numerical  method  of  solution  section. 

It  Is  important  to  note  that  this  approximation  retains  all  of  the  terms 
which  would  arise  in  a  second-order  boundary-layer  theory. 
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It  can  be  shown  that  with  the  parabolic  approximation,  the  stream- 
function  Equation  (2.5)  reduces  to  the  boundary-layer  form  and  is  given 


The  steady-state  boundary-layer  equations  in  optimal  coordinates  are 
obtained  from  the  full  Navier-Stokes  Equations  (2.5)  and  (2.12)  by  taking 
the  limit  of  the  equations  as  Reynolds  number  Re  -*■  <»,  in  inner  boundary- 
layer  variables.  By  making  the  boundary-layer  approximation  in  Equations 
(2.5)  and  (2.12)  it  can  be  shown  that  many  terms  drop  out  and  in  addition 
the  remaining  coefficients  h,  H^,  y^,  and  can  be  evaluated  at  the  body 
surface,  i.e.  at  n  *  0.  These  coefficients  can  therefore  be  written  as 
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h  -  h(£) 


J  J  (5) 

r  r 


and 


Yr  Yr«) 


H  -  t  (5,0) 

r 


(2.31a-c) 


With  these  approximations  the  streamf unction  Equation  (2.5)  remains  the 
same,  i.e. 


f  »  g 
nn  6 


(2.32) 


and  the  vorticity  Equation  (2.12)  reduces  to 


8nn  +  ~fr  <f+5V8n  '  £  (1  +  5  iK1*  '  S  £  *  0  (2'33> 

The  inverse  scale  factor  H  in  the  case  of  boundary-layer  can  be  written  as 


H(5,0) 


\  ,  t2(5,0) 
h2  h2(C,0) 


(2.34) 


The  expression  for  the  scale  factor  h  is  given  once  the  geometry 
is  prescribed  and  will  be  given  later  in  the  respective  example  problem. 

Now  substituting  f  -  g  into  the  vorticity  Equation  (2.33)  and 
integrating  the  resulting  equation  with  respect  to  n  and  evaluating  the 
constant  of  integration  by  making  use  of  the  boundary  condition  on  f 
and  its  derivatives  as  n  -*■  °»,  the  vorticity  equation  can  be  written  as 


A 
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J.  J  O  J  ,  H,  „ 

f  +  —  (f+5f  )f  -  —  f 2  -  —  4  f2 

nnn  Yr  5  nn  Yr  n  Yr  2  h  n 


J  J  _  H, 

_  r  —  f  f  _ _ -(!+-§-  _£) 

4  Y_  Cn  Y.  U  2  H  J 


(2.35) 


We  could  work  with  these  variables,  however  It  is  convenient  to  instead 
develop  variables  similar  to  those  of  Gortler.  In  order  to  do  this  we 
define  V,  a  normal  velocity  like  variable,  as 


V  =  -  f  -  £f , 


(2.36) 


Making  suitable  rearrangements  in  the  Equation  (2.35)  and  replacing  f 


by  F  we  get  the  final  form  of  the  boundary-layer  equations  as 


nn 


and 


J  J  .  J 

F - —  VF  +  —  8.  (1  -  F  J  -  £  —  FF_ 

—  Y  n  Y  i  Y  £ 

r  fr  r 


-  F  -  £F 


n  -'C 

where  the  pressure  gradient  parameter  6^  is  given  by 

CH 

Si  -  d  +  2H  > 


(2.37) 


(2.38) 


(2.39) 


In  transformed  variables,  a  displacement  thickness  like  variable  6 
can  be  defined  as 


<5e*  /  (1  -  F)dn  (2.40) 

o 

The  necessary  boundary  conditions  to  solve  the  governing  differential 
Equations  ((2.38)  and  (2.39))  are  as  follows.  The  no-slip  and  zero-in¬ 
jection  conditions  at  n  *  0  give 


21 


F(£,0)  =  0 


and  (2.41a,b) 

V(£,0)  -  0 

The  asymptotic  boundary  condition  in  the  transverse  direction  gives 

F(£,n)  ^1  as  n  -*•  ®  (2.42) 

The  flow  is  symmetric  about  the  line  £  *  0,  therefore  the  symmetry  con¬ 
dition  can  be  written  as 

F?  *  0  at  £  -  0  .  (2.43) 

For  the  bodies  we  will  be  considering,  at  downstream  infinity  J^/y^  *►  1 

and  the  pressure  gradient  parameter  8  ^  approaches  the  value  of  the  Falkner- 

Skan  pressure  gradient  parameter  8  and  the  equation  reduces  to  the  Falkner- 

Skan  equation  which  can  be  rewritten  as 

F  -  VF  +  8(1  -  F2)  =  0  (2.44) 

nn  n 

and 

V  -  -  F  (2.45) 

n 

It  is  interesting  to  observe  that  the  boundary-layer  equations  in  the 

present  coordinates  show  self-similar  behavior  and  similarity  solutions  are 

automatically  recovered  whenever  similarity  exists.  It  is  also  important 

to  remark  that  the  present  form  of  the  boundary-layer  equations  is  similar 

** 

to  those  when  written  in  Gortler  variables. 

This  completes  the  general  analysis  for  determination  of  the  boundary- 
layer  equations  in  optimal  coordinates  except  for  the  missing  boundary 
conditions  t(£,0).  In  the  following  we  will  develop  a  scheme  to  determine 
the  surface  boundary  condition  on  t(£,n). 
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5 .  The  Surface  Boundary  Condition  for  the  Boundary -Layer  Equations 


In  the  course  of  our  previous  discussion  it  was  mentioned  in  Section 
II  that  the  coordinate  equation  is  coupled  to  the  flow  equation  through 
the  surface  boundary  condition  on  t(£,n).  The  optimal  coordinate  con¬ 
dition  states  that  the  transformed  displacement  thickness  C  should  be 
constant.  We  will  make  use  of  this  condition  to  determine  the  missing 
boundary  condition  on  t(£,n).  Since  there  is  no  explicit  analytical 
relation  relating  the  surface  boundary  condition  on  t(£,0)  and  the 
displacement  thickness,  C  in  optimal  coordinates,  we  must  determine  the 
value  of  t(£,ri)  on  the  surface  iteratively.  The  iterative  relation  which 
yields  the  value  of  t(£,0)  in  optimal  coordinates  is  determined  in  the 
following  manner  along  with  some  other  interesting  properties  of  the 
equations. 

Let  us  assume  £  and  n  are  the  coordinates  that  are  optimal  and  hence 
are  to  be  determined.  The  coordinates  £  and  n  are  assumed  not  to  be 
optimal  but  some  suitable  coordinate  system  which  is  used  as  an  initxal 
guess  for  the  optimal  coordinate  system.  Let  us  also  assume  that 

f  -  A(£)f(£,n)  ,  (2.46) 

n  -  A(£)n  (2.47) 

and 

£  -  B(£)  ,  (2.48 

where  A  and  B  are  arbitrary  functions  of  £ . 
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Then  it  follows  that 


24 


and 


(2.55) 


Combining  the  relations  given  by  Equations  (2.52)  and  (2.55)  we  get 


A2f  (1+5fl)  =  —  —  A^B ' 

T  A  l  \ 


which  can  be  simplified  to  yield 


(2.56a) 


=  1  =  B(C)  ‘  (2.56b) 

A 

The  above  relation  is  already  satisfied  by  Equation  (2.53).  There¬ 
fore  we  can  ignore  one  of  the  two  relations  given  by  Equations  (2.52) 
and  (2.55).  Hereafter  we  will  neglect  the  relation  (2.55). 

Equation  (2.56b)  gives 


5B' 


B(5) 


or  rearranging  we  obtain 


(2.57a) 


Integrating  both  sides  the  Equation  (2.57a)  we  get 

B  *  K5A  (2.57b) 

where  K  is  a  constant  of  integration  which  will  be  determined  later. 


25 


Using  the  relation  (2.57b)  and  (2.53)  we  can  write 
£A'  B ' 

(1+a)=Ia  *  (2-58) 

Substituting  the  relation  (2.58)  into  the  Equation  (2.52)  we  get 


Jr  _  AB’  Jr 
-  K  Y_ 


(2.59) 


and  from  the  relation  (2.57b)  and  (2.52)  we  can  write 
l  -  K?A 


(2.60) 


Now  substitution  of  J^/y^  from  Equation  (2.59)  into  Equation  (2.54)  yields 


a2  Jr  .  AB'  Jr  - 
Yr  6i  "  K  Yr  61 


Simplifying  the  relation  (2.61)  we  obtain 


(2.61) 


0. 


AK 


3i  "  B'  Bi 
where  Bi  is  given  as 

%  ■ 


(2.62) 


(2.63) 


which  can  be  rewritten  as 


%'■ 


if!?** 

(H52)  2 


(2.64) 


Using  the  relations  (2.58)  and  (2.64)  in  the  Equation  (2.62)  we  can  write 
8*  as 


8 


i 


(C2H)g 

(52H) 


(f) 


(1  +541) 


(52H) 

(52H) 


I  I  A 
2  «A)» 


(2.65) 
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Also  we  have  from  Equation  (2.60) 


C  is  related  to  the  displacement  function  in  optimal  coordinates. 

In  general  in  the  (£,n)  coordinate  system  (which  may  not  be  optimal) 
the  streamf unction  f(£,n)  far  from  the  body  is  given  as 

f  *  n  -  US)  (2.69) 
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with  the  help  of  relations  (2.45)  and  (2.46)  the  Equation  (2.69)  can  be 
written  as 


f (f.n) 


^  n 


Kg) 

A(C) 


as 


n  -*■  » 


(2.70) 


Comparing  the  two  expressions  for  f (C,n)  given  by  Equations  (2.68)  and 


(2.70)  we  get 

c  ■  sfj-  ■ 1  <2-71> 

Also  we  can  write  from  Equation  (2.71) 

A(C)  =  Aiil  .  (2.72) 

Substituting  the  expression  of  A(?)  from  Equation  (2.71)  and  its  derivative 


in  the  relations  (2.59)  through  (2.61)  we  get 


Kg)(gKg))* 


(2.73) 


g 


.and 


H.HO 

i 


(2.74) 


H 


(gKg))' 


(2.75) 


Also  we  can  write  in  the  present  coordinate  system  for  the  boundary- 

layer  flow  from  Equation  (2.25d)  as 


Yr  t(t+Ct^) 


(2.76) 


Similarly  we  can  write  an  analogous  expression  for  barred  quantities  as 


3 

r 


1 

t(t+Ctc) 


(2.77) 
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From  the  coordinate  relation  (2.23)  we  can  write 

p  =  It  -  It  (2.78) 

Substituting  \  and  5  from  the  Equation  (2.78)  into  the  Equation  (2.74) 
we  get 


-  «  K  (2.79) 

t  £ 

Using  the  relations  (2.76)  and  (2.77)  in  Equation  (2.73)  we  can  show  that 


_  *(S)(g*(5))» 


tG+ll) 


-2  t(t+5t?) 


(2.80a) 


or 


K(— )  =  (-)  (^)  1 


£  (£t)^  £  £ 


(2.80b) 


Now  using  the  relation  (2.78)  we  can  easily  write,  after  rearranging  the 
equation  (2.80b)  as 

K  =  d(5p)  I  d£ 
d5  £  dC 


or 


(2.80c,d) 


d|  d(gp)  1  d£ 
d5  dS  £  d? 

which  can  be  simplified  to  give 


di.ll  d(S£) 
d5  K  -  d5 

In  light  of  the  relation  for  £  given  in  Equation  (2.74), 
C2,8Qe)  is  true  if  and  only  if  K  ■  1, 


(2.80e) 
the  relation 
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Therefore  the  appropriate  relation  and  the  one  which  can  be  used 
is 


where  t  is  now  the  needed  boundary  condition  on  the  surface  for  the  co¬ 
ordinate  function  t(£,n)  and  l  is  the  displacement  thickness  function  in 
optimal  coordinates.  The  quantities  t  and  £(£)  are  the  corresponding 
variables  in  the  guessed  (£,n)  coordinate  system. 

This  completes  the  brief  outline  of  the  development  of  the  governing 
differential  equations  together  with  the  required  boundary  conditions  for 
the  two  dimensional  symmetric  incompressible  flow  in  optimal  coordinates. 
The  details  of  the  numerical  analysis  necessary  to  solve  the  given  set  of 
equations  developed  in  this  section  are  elaborated  on  in  the  the  follow¬ 
ing  section. 
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III.  NUMERICAL  METHODS  OF  SOLUTION 


1 .  Boundary-Layer  Solution  Method 

A  marching  procedure  is  made  use  of  in  determining  the  solution  of 

the  boundary-layer  equations  in  optimal  coordinates.  The  only  difference 

between  the  present  method  and  the  boundary-layer  forward  marching  approach 
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used  by  numerous  authors  (see  Blottner  and  Flugge  Lotz  for  example)  is 

that  an  additional  iterative  loop  has  been  introduced  here  in  order  to 

satisfy  the  condition  of  optimal  coordinates.  The  coefficients 

a  ,  8  ,  Y  and  J  which  are  functions  of  t,  are  evaluated  at  the  surface 
r  r  r  r 

if  they  are  present  in  the  boundary-layer  equation.  This  fact  represents 
an  advantage  of  setting  the  iterative  scheme  for  the  optimal  coordinate 
condition  without  solving  the  coordinate  equation  in  each  iteration. 

This  implies  that  the  surface  boundary  condition  t(5,0)  on  the  coordinate 
Equation  (2.25)  can  be  computed  once  and  for  all  by  iterating  the  solution 
vectors  of  the  boundary-layer  equations  such  that  the  transformed  displace¬ 
ment  thickness  function  C  is  constant  within  a  prespecified  limit  of 
accuracy.  To  compute  t(5,0)  in  each  iteration  the  relation  (2.81)  is 
used  to  update  the  surface  boundary  condition  t (5,0). 

In  order  to  solve  the  boundary-layer  equations,  a  coordinate  trans- 
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fromation  similar  to  that  used  by  VandeVooren  and  Dijkstra  ,  Davis 
and  Davis,  U.  Ghia  and  K.  Ghia^  is  used  in  the  5(streamwise)  direction 
to  achieve  proper  resolution  of  the  flow  field  in  the  streamwise  direction. 
The  transformation  which  gives  the  appropriate  flow  field  resolution  is 


given  by 

S  -  1  -  A/5  log(l  +  5/A) 


(3.1) 


where  A  is  a  constant,  the  value  of  which  is  chosen  such  that  the  vorti- 
city,  g,  at  the  wall  drops  to  one  half  the  stagnation  point  value  at 
about  S  =  0.5,  where  S  is  the  transformed  coordinate.  The  value  of  A 
which  gives  such  a  variation  of  vorticity  g(C.n)  at  the  wall  is  found  to 
be  equal  to  1  for  the  boundary-layer  analysis.  Apart  from  the  advantages 
mentioned  above,  the  transformation  relation  (3.1)  transforms  the  semi¬ 
infinite  domain  of  integration  in  the  streamwise  (£)  direction  into  a 
finite  domain  in  the  new  independent  variable  t  such  that  S  varies  between 
0  to  1  as  5  varies  from  0  to  ®. 

It  is  mentioned  by  Davis^  that  the  vorticity  g  dies  out  exponentially 
to  zero  beyond  a  value  of  n  =  5.0  for  all  Reynolds  numbers  in  the  case  of 
laminar  unseparated  flow.  The  boundary-layer  equation  is  therefore  inte¬ 
grated  up  to  n  =  6.0  with  an  uniform  step  size  of  An. 

In  order  to  solve  the  boundary-layer  equations,  the  momentum  equation 
which  is  coupled  to  the  normal  velocity  equation  is  linearized  with  res¬ 
pect  to  the  previous  step  value  and  rewritten  in  new  independent  vari¬ 
ables  as 


F 

nn 


-  v,0V(\  Vg(1  -  V>  -  ^  V  at 


Then  the  normal  velocity  equation  reads 
V  +  F  +  ZI  ||  -  0 

n  s  dc 


(3.2) 


(3.3) 


where  the  subscript  G  stands  for  the  guessed  value  of  the  dependent  vari¬ 
ables  from  the  previous  station  or  iteration. 
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An  implicit  finite  difference  scheme  similar  to  that  given  by 
22 

Blottner  and  Flugge-Lotz  is  used  to  solve  the  parabolic  momentum 

Equation  (3.2).  Once  the  momentum  equation  is  solved  the  continuity 

equation  is  integrated  with  respect  to  n  once  to  obtain  the  updated 

value  of  V  to  be  used  in  the  subsequent  iteration  of  the  momentum 

equation.  In  the  numerical  scheme,  the  F  and  F  derivative  terms 

n  nn 

are  written  in  second  order  accurate  central  difference  form  and  the  S 

derivative  is  written  as  a  two  point  backward  difference.  Application 

of  such  a  scheme  to  the  momentum  Equation  (3.2)  results  in  a  tria- 

diagonal  system  of  equations  in  the  p  direction  at  each  grid  line  given 

24 

by  S  =  constant.  The  Thomas  algorithm  is  employed  to  invert  the  set  of 
simultaneous  algebraic  equations.  The  momentum  and  the  continuity 
equations  are  solved  simultaneously.  The  solution  vector  is  iterated  on 
due  to  the  nonlinearities  until  the  values  of  the  solution  obtained  in 
two  successive  iterations  do  not  differ  by  a  prespecified  accuracy  limit. 
Repeating  these  steps  at  each  grid  line  at  S  =  constant,  the  solution 
procedure  is  marched  from  S  *  0  to  S  «  1  to  cover  the  entire  semi-infinite 
domain  in  the  streamwise  direction.  In  order  to  initiate  the  marching 
process  we  need  the  initial  profile  of  the  solution  vectors  F  and  V  at 
5*0  and  also  the  value  of  the  displacement  function  1  in  optimal  coordi¬ 
nates  in  order  to  calculate  t(5,0)  along  the  surface  of  the  body.  The 
displacement  thickness  in  optimal  coordinates  and  the  starting  values 
of  the  solution  vectors  F  and  V  at  the  stagnation  point  are  determined 
in  the  following  manner. 
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It  was  already  mentioned  in  the  Section  II  on  the  asymptotic  boundary 
conditions  for  the  coordinate  Equation  (2.25)  that  the  optimal  coordinates 
become  conformal  at  downstream  infinity.  This  Implies  that  the  displace¬ 
ment  thickness  function  £(£)  at  downstream  infinity  is  the  displacement 
thickness  function  l  in  optimal  coordinates.  Therefore  in  order  to  de¬ 
termine  the  displacement  thickness  function  7  the  boundary  layer  Equations 
(3.2)  and  (3.3)  are  solved  at  downstream  infinity.  As  5  ^  the  boundary- 
layer  equations  yield  the  similarity  form  as  given  below. 

F  -  VF  +  0 (1  -  F2)  -  0  (3.4) 

Tin  n 

and 

Vn  +  F  -  0  (3.5) 

where  8  assumes  a  value  equal  to  zero  for  the  bluned  plate  case  and 

8  *  25/v  for  the  blunted  wedge  problem. 

The  above  set  of  Equations  ( (3 .4) — (3 . 5) )  can  easily  be  solved  with 

the  boundary  condition  prescribed  by  Equations  (2.41a,b)  and  (2.42) . 

Now  writing  the  derivatives  P  and  F  in  central  difference  form  a  set 

nn  n 

of  triadiagonal  system  of  equations  are  obtained  and  these  can  easily  be 
inverted  by  applying  Thomas  algorithm.  Once  the  solution  is  obtained, 
the  displacement  thickness  I  is  calculated  using  the  Equation  (2.40). 
Finally,  the  initial  profiles  of  P  and  7  are  obtained  by  solving  the 
boundary-layer  equations  at  5  ■  0,  Again  we  note  from  Equations  (3.2) 
and  (3,3)  at  (  ■  0,  that  the  boundary-layer  equations  assume  a  similarity 
form  because  of  the  symmetric  condition  F^  -  0  along  the  the  stagnation 
line.  The  equations  at  the  stagnation  line,  i.e.  at  ?  ■  0,  reduce  to 
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(3.6) 


O  U  « 

F - -  VF  +  —  (1  -  F  )  =  0 

nn  rr  n  Yr 

and 

V  +  F  *  0  (3.7) 

n 

with  the  boundary  conditions 
F  *  0 

at  n  =*  0  (3.8a,b) 

V  -  0 

and 

F  ^  1  as  n  +  ®  .  (3.8c) 

The  Equations  (3.6)  and  (3.7)  can  easily  be  solved  in  exactly  the  same 
manner  as  in  the  case  of  downstream  infinity.  But  in  this  case  the 
solution  of  the  boundary-layer  Equations  (3.6)  and  (3.7)  has  to  be  iterated 
on  to  satisfy  the  condition  of  optimal  coordinates,  in  addition  to  the 
iteration  required  for  the  nonlinearity  of  Equation  (3.6).  Once  the  init¬ 
ial  profile  is  obtained  it  is  now  possible  to  compute  the  solution  of  the 
Equations  (3.2)  and  (3.3)  in  the  entire  flow  field  by  using  a  forward 
marching  process.  After  the  determination  of  t(£,0)  at  all  of  the  grid 
points  on  the  surface,  the  coordinate  Equation  (2.25)  can  be  solved. 

The  details  of  the  numerical  analysis  for  the  coordinate  Equation  (2.25) 
is  discussed  in  the  following. 

2.  Optimal  Coordinate  Solution  Method 

The  coordinate  equation  as  given  by  (2.25)  is  elliptic  in  nature 
and  requires  four  boundary  conditions  of  which  three  are  already  specified 
by  the  Equations  (2.27)  through  (2.28a,b).  The  fourth  boundary  condition 
is  determined  from  the  flow  equations  as  mentioned  in  the  preceding  section. 
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Because  of  the  elliptic  nature  of  the  equation  an  alternating  direction 

25 

implicit  (ADI)  method,  first  developed  by  Douglas  ,  is  made  use  of  to 

obtain  the  coordinate  equation  solution.  In  order  to  apply  the  ADI 

scheme  a  fictitious  unsteady  term  3t/3v  is  added  to  the  right  hand  side 

of  the  Equation  (2.25).  The  resulting  equation  is  marched  in  time  and 

alternated  with  respect  to  the  independent  variables  C  and  nnat  each  half 

time  step.  The  steady  state  solution  is  recovered  in  a  two  step  process 

to  obtain  the  asymptotic  converged  solution  of  the  unsteady  differential 

equation.  Recognizing  the  fact  that  the  transient  solution  is  of  no 

importance  for  the  present  analysis,  certain  measures  are  taken  in  order 

to  achieve  faster  convergence.  For  example,  the  coefficient  of  the  term 

3t/3v  is  set  equal  to  different  constants  for  different  problem  parameters 

to  improve  the  convergence  rate  of  the  ADI  scheme. 

A  transformation  of  the  independent  variables  is  made  in  order  to 

confine  the  domain  of  integration  to  a  finite  region  besides  providing 

the  proper  resolution  in  the  n  and  5  directions  and  preserving  the 
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uniform  step  size  accuracy  in  the  computational  domain.  The  trans¬ 
formation  relations  which  serve  these  purposes  are  given  in  the  following 
for  the  n  variable  and  in  (3.1)  for  the  5  variable.  The  variable  n  is 
replaced  by 


In  Equation  (3.1)  the  constant  A  as  mentioned  previously  is  chosen  equal 
to  1.  The  transformation  (3.1)  together  with  (3.9)  transform  the  quarter 
infinite  domain  in  (C,n)  coordinates  to  a  unit  square  in  the  computational 


domain  of  (S,N)  coordinates.  Using  these  transformation  relations  the 
coordinate  Equation  (2.25)  is  rewritten  in  new  independent  variables 
(S,N)  as 

2  ,2„  23 


t-  (M)  +  r  a  t  (— )  (— )  +  — 

T*N  W  TJ  .  2  y  CSN  W  W  y 
dp  'r  'r 


d2N 


(1  dS  +d^S  +^r  dS  2  =  3t 

kE,  d?  dc2  S  Yr  SS  '■dr  sv 


(3.10) 


The  subscripts  S,N  in  the  Equation  (3.10)  represent  the  partial  differ¬ 
entiation  of  the  dependent  variable  with  respect  to  the  independent 
variables  appearing. 

The  cross  derivative  term  t  in  the  coordinate  equations  needs  a 

O  IN 

special  treatment  because  it  plays  an  important  role  in  the  numerical 
solution  of  the  coordinate  equation.  The  cross  derivative  term  appears 
in  the  equation  because  of  the  non-orthogonality  of  the  (£,n)  coordinate 
system.  The  coordinate  system  depends  on  the  displacement  effect.  The 
magnitude  of  this  cross  derivative  term  is  also  an  indication  of  the 
degree  to  which  the  optimal  coordinates  are  non-orthogonal,  i.e.  vary  from 
conformal  coordinates  in  this  case.  One  must  also  be  careful  in  formu¬ 
lating  the  numerical  scheme  for  the  cross . derivative  since  it  is  well 

known  that  the  cross  derivative  if  improperly  handled  can  cause  numerical 
27 

instability.  In  many  situations  this  term  is  small  and  can  be  treated 
28 

explicitly  in  an  otherwise  implicit  scheme.  In  the  present  analysis, 
in  order  to  be  consistent  with  the  numerical  scheme,  an  implicit  treatment 
of  the  cross  derivative  has  been  developed  using  a  three  level  finite 


difference  scheme  similar  to  that  followed  by  Mitchell 
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This  type  of 
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the  treatment  of  the  cross-derivative  term  allows  one  to  break  the  single 
term  t_^  into  combinations  of  unidirectional  derivatives  with  respect  to 
the  independent  variables  £  and  n.  The  cross  derivative  t  is  written 
using  the  three  level  grid  points  as  shown  in  Figure  2.  Then  the  cross¬ 
derivative  at  the  grid  point  m,n  can  be  expressed  in  second  order  accurate 
difference  form  as 


t  -  I-DF  ~  t,  +  op  ££  t 

5n  AS  TWi  +  DF  AS  SfNt  . 

|m,n  |m-l,n 


,  1  _1_  AS 

2AS  CN,  ‘  2AS  \  .  ~  DF  AN  CSS , 

I  m-1 ,  n  |m,n 


2DF  ,  (DF  "  AS  . 

AN  S,m>n  2  An'sS,^ 


(DF  -  i)  (DF  +  h  ... 

+ - c  + _ L.A2.* 

AN  S.  2  AN  tSS, 

|m,n-l  |m,n+i 


(DF  +  i) 


(3.11) 


The  other  terms  in  Equation  (3.11)  are  written  as  three  point  central 
differences.  It  can  be  shown  that  in  order  to  maintain  diagonal  dominance 
in  the  inversion  scheme,  DF  should  be  taken  equal  to  0,  1/4  or  -  1/4 
according  to  whether  the  coefficient  of  the  t^  term  in  Equation  (3.10) 

^ 8  zero»  positive  or  negative.  Keeping  track  of  this  property  provides 
stability  to  the  system.  The  ADI  numerical  method  of  solution  using  this 
difference  scheme  is  described  in  the  following. 
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During  the  first  half  time  step,  denoted  by  the  star  (*)  step,  the 

ADI  algorithm  for  the  coordinate  Equation  (3.10)  is  advanced  from  time 
n  *  n  fay 

v  to  v  =  v  +  ~2~  by  following  relation: 


8  8 
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(3.12) 


where  Av  is  the  time  step. 
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The  starred  (*)  quantities  in  Equation  (3.12)  are  considered  to  be 
unknown  whereas  the  dependent  variables  with  superscript  n  are  assumed 
to  be  known  from  the  preceding  half  time  step  (n) .  The  Equation  (3.12) 
is  solved  along  lines  of  constant  S  or  5.  Writing  the  t^  and  t^  terms 
in  central  difference  form,  the  Equation  (3.12)  can  be  shown  to  reduce 
to  a  set  of  simultaneous  equations  in  triadiagonal  form.  These  algebraic 
equations  can  be  inverted  by  applying  the  Thomas  algorithm. 


During  the  second  half  time  step.  Equation  -(3.12)  is  marched  forward 
n.-t'l  ^  Av 

to  a  new  time  level  v  =  v  +  y-  using  the  relation  given  below: 
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(3.13) 
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where  the  starred  (*)  quantities  are  known  from  the  previous  half  time 
step  calculation  and  the  dependent  variables  with  superscript  (n+l)  are 
considered  to  be  unknown  in  the  Equation  (3.13).  Writing  the  derivatives 
tgs  and  tg  in  second  order  accurate  central  difference  form,  the  differential 
Equation  (3.13)  can  be  solved  in  an  identical  manner  to  the  star  (*)  step.  The 
solution  procedure  for  the  coordinate  equation  alternates  between  the 
constant  £  and  constant  n  directions  in  one  complete  time  step.  The 
steady  state  solution  is  achieved  as  an  asymptotic  converged  solution 
to  the  unsteady  Equation  (3.10).  Since  the  initial  guessed  coordinates 
assumed  for  the  iteration  procedure  are  assumed  to  be  conformal,  the 
initial  conditions  for  t(£,n)  at  each  grid  point  is  taken  to  be  1. 

Next  we  discuss  the  numerical  solution  of  the  Navier-Stokes  equations 
in  optimal  coordinates.  However,  rather  than  solving  the  complete  equa¬ 
tions  we  solve  a  parabolized  version.  The  solution  of  the  complete  equa¬ 
tions  should  follow  in  a  similar  manner. 

3.  Solution  of  Parabolized  Navier-Stokes  Equations 

The  parabolized  vorticity  and  streamf unction  equations  unlike  the 
boundary-layer  equations,  are  coupled  with  the  coordinate  equation  because 
the  coeff icients  af,  6^,  Yr  and  are  no  longer  functions  of  £  only  but 
are  functions  of  n  as  well.  This  fact  calls  for  the  simultaneous  treatment 
of  the  flow  Equations  ((2.29)and  (2.30))  and  the  coordinate  Equation  (3.10). 
Therefore,  the  iterative  scheme  described  in  the  case  of  boundary-layer 
analysis  is  modified  slightly  to  accommodate  the  simultaneous  treatment 

of  the  flow  equations-  and  the  coordinate  equation.  Details  of  the  scheme  ] 

< 

are  presented  in  the  following, 

t 

* 

i 
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In  order  to  solve  the  vorticity  and  streamf unction  equations  more 
efficiently  transformations  of  the  independent  variables  are  made.  The 
transformation  relations  are  the  same  as  given  in  case  of  the  coordinate 
equation  and  are  given  by  the  Equations  (3.1)  and  (3.9).  The  only  differ¬ 
ence  in  the  use  of  the  relation  (3.1)  for  the  case  of  parabolized  Navier- 
Stokes  analysis  is  in  the  value  of  the  constant  A,  which,  in  this  case 
is  chosen  to  be  4.0  +  0.4  q^,  where  q^  =  q^  /Re  .  Re  is  the  Reynolds 
number  based  on  the  nose  radius  of  curvature,  a  and  q^  is  the  value  of 
q  at  the  wall.  The  dependence  of  the  constant  A  on  Re,  in  the  stream- 
wise  stretching  is  introduced  in  order  to  achieve  the  desired  flow  field 
resolution  in  the  streamwise  direction  for  different  Reynolds  numbers. 

The  transformation  relation  (3.9)  in  the  n  direction  has  a  viscous  type 
stretching  built  into  it,  such  that  it  always  places  the  viscous  layer 
in  approximately  the  lower  half  of  the  unit  square  region  of  integration 
while  the  inviscid  flow  is  contained  in  the  upper  half.  This  type  of  flow 
field  resolution  in  the  n  direction  is  appropriate  in  light  of  the  fact 

that  for  all  unseparated  flow  cases,  as  was  mentioned  in  the  case  of 

% 

boundary-layer  analysis,  the  viscous  effects  are  confined  within  the 
region  n  <  5  or  in  the  transformed  variable  N  <  0.5.  The  transformation 
relation  (3.9)  therefore  allows  one  to  obtain  proper  grid  spacing  in  the 
physical  flow  field  near  the  wall,  where  the  gradients  of  the  flow  vari¬ 
ables  are  large  as  compared  to  their  asymptotic  values  far  away  from  the 
wall.  Various  other  advantages  of  the  transformation  relations  (3.1) 
and  (3.9)  are  elaborated  on  in  the  case  of  boundary-layer  analysis  and 
will  not  be  repeated  here. 
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Because  of  the  boundary  condition  f+nrCasn+®,  a  new  stream- 
function  h  is  defined  such  that 

f  =  h  +  n  (3.15) 
The  steady  state  streamfunction-vorticity  equations  in  the  transformed 
variables  can  now  be  written  as 
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a  S,H  B  N 

a  +  +  I2^iP  -  2  h’  0 


(3.17) 


The  corresponding  boundary  conditions  are  given  in  the  following.  The 
no  slip  and  zero  injection  conditions  on  the  surface  give 


h  -  0 


and 


at  n  •  0 


(3.18a,b) 
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The  infinity  boundary  conditions  are  given  by 

h  v  0 

n 

and  as  n  -*■  «  (3.19a,b) 

g(5,n)  0 

At  C  *  0,  symmetry  is  made  use  of  to  compute  the  vorticity  and  streamf unction 
along  the  stagnation  plane. 

The  asymptotic  boundary  condition  as  £  -*•  “  is  given  by  Equation  (3.16) 
along  with  (3.17)  evaluated  as  £  ->•<».  Equation  (3.16)  remains  the  same 
whereas  (3.17)  reduces  to 

(Nn)2gNN  +  <1  +  S  +  1)Nn%  +  <!  "  26)  (1  +  hNNri)g  (3.20) 

where  6  is  the  Falkner-Skan  pressure  gradient  parameter  and  is  a  known 
quantity  for  a  given  problem. 

It  is  observed  from  Equations  (3.16)  and  (3.17)  that  the  streamf unction 
and  the  vorticity  equations  are  second-order  partial  differential  equations 
with  respect  to  n.  Therefore  two  boundary  conditions  for  each  of  the 
dependent  variables  h  and  g  are  needed  in  order  to  integrate  the  equations. 
There  are  three  boundary  conditions  on  h  and  one  boundary  condition  on  g. 
These  type  of  boundary  conditions  are  not  the  most  convenient  ones  for  the 
Thomas  algorithm.  However,  the  solution  to  the  flow  equations  can  be 
obtained  by  a  superposition  method  as  described  by  Davis,  U.  Ghia  and 
K.  Ghia, ^  by  determining  a  coe-ficient  such  that 

g(£,n)  -  c i8h(  £,n  5  +  gp(£,n)  (3.21) 
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with  given  by 


(3.22) 


where  g  and  g  are  the  homogeneous  and  particular  solutions  respectively 
h  p 

of  the  vorticity  Equation  (3.17)  and  h  and  h  are  the  particular  solu- 

P1  p2 

tions  of  the  streamf unction  Equation  (3.16)  corresponding  to  the  forcing 

function  g  (£,n)  and  g  (5,n). 

H  p 

The  governing  flow  equations  are  therefore  integrated  along  the 


curves  of  constant  £  for  all  n,  using  the  implicit  finite  difference 

**  22 

method  of  Blottner  and  Flugge-Lotz  as  described  in  the  case  of  boundary- 
layer  analysis.  Once  the  required  convergence  is  achieved  on  the  solution 
vectors,  the  t(5,0)  value  is  updated  using  the  iterative  relation  (2.81). 
With  this  surface  boundary  condition  distribution  on  t,  the  coordinate 


equation  is  Integrated  employing  the  two  step  time  dependent  ADI  relaxation 

method,  as  described  previously  in  boundary-layer  analysis.  The  converged 

solution  is  then  stored  to  calculate  the  new  values  of  J  ,  a  ,  8  ,  and  H 

r  r  r  r 

to  be  used  subsequently  in  the  flow  equations.  Before  starting  a  new  iter¬ 
ation  in  the  flow  equations  after  this  step,  the  edge  quantity  h(5,°°)  is 
checked  for  satisfaction  of  the  condition  of  optimal  coordinates,  i.e. 

1(0  -*■  C,  a  constant.  If  this  condition  is  not  satisfied  within  a  pres¬ 


cribed  accuracy  limit,  then  the  entire  cycle  is  repeated  to  calculate  a 
new  set  of  values  for  the  solution  vectors  g,  f  and  t. 
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This  completes  the  method  of  solving  the  flow  equations  in  optimal 
coordinates  including  the  generation  of  the  optimal  coordinates.  The 
application  of  the  analysis  to  specific  problems  is  discussed  in  the 
following  section. 
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IV.  APPLICATION  TO  SYMMETRIC  TWO-DIMENSIONAL  FLOW  PROBLEMS 


In  order  to  demonstrate  the  application  of  the  principles  described 
in  the  previous  sections  we  will  generate  the  solutions  for  two- 
dimensional  incompressible  laminar  symmetric  flow.  Two  geometries  are 
considered  as  examples;  the  flow  past  a  semi-infinite  wedge,  and  the 
flow  past  a  semi- inf inite  thick  flat-plate,  see  Figures  3a, b  and  4  and  5 
respectively.  Boundary-layer  as  well  as  the  parabolized  vorticity- 
streamfunction  form  of  the  flow  equations  are  used  for  the  present 
analysis.  Details  of  the  analysis  are  presented  in  the  following. 


1.  The  Blunted  Wedge  Example 

The  two  dimensional  incompressible  flow  past  a  semi-infinite  wedge 
as  shown  in  the  Cartesian  plane  in  Figure  3a  is  considered  here.  The 
axis  of  the  wedge  is  oriented  in  the  direction  such  that  the  flow  is 
symmetric  about  the  wedge  axis.  The  body  surface  in  the  physical  plane 
as  shown  in  Figure  3a  is  mapped  to  the  line  q  ■  0  and  the  corresponding 
flow  in  the  physical  (x,y)  plane  is  mapped  into  a  stagnation  point 
type  flow  which  is  the  upper  half  of  the  ^  plane  (see  Figure  3b).  This 
type  of  mapping  is  achieved  by  applying  a  Schwarz-Christoffel  trans- 


(4.1) 


where  5  ■  q  +  ip.  The  scale  factor  of  the  transformation  is  given 


(4.2) 
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The  factor  v/U^  appearing  on  the  right  hand  side  of  Equation  (4.1), 
non-dimensionalizes  the  coordinates  q  and  p.  The  n,  appearing  in  the 
exponent  of  the  transformation  relation  (4.1)  is  related  to  the  semi¬ 
wedge  angle  6  as 


6 

n-6 


(4.3) 


Recognizing  that  the  point  (yQ,0)  transforms  to  the  point  (^,0), 

(see  Figures  3a  and  3b),  the  nose  radius  of  curvature  a  can  be  written, 
using  Newton's  formula,  as 


,1+iu  1+n  v 

*  *  W  s,  r 


(4.4) 


where  q  is  equal  to  q^  at  the  wall.  Defining  the  Reynolds  number  based 

on  the  nose  radius  of  curvature  such  that  Re  ■  U  a/v,  it  can  be  shown 

r 

from  the  relation  (4.4)  that  the  location  of  the  wedge  surface  in  the 
transformed  plane  is  given  by 


1+n 
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(4.5) 


It  is  observed  from  the  equation  (4.5)  that  keeping  the  Reynolds 
number  Re  constant,  the  effects  of  body  shape  on  the  flow  can  be  studied 
by  changing  n,  whereas  keeping  n  fixed,  the  effects  of  Reynolds  number 
on  the  flow  can  be  investigated  by  varying  Re. 


m 
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It  will  be  shown  later  that  the  conformal  transformation  used 
here  to  transform  the  physical  flow  field  into  a  stagnation  point  flow 
in  the  conformal  plane,  allows  one  to  obtain  similarity  type  variables 
and  similarity  solutions  are  recovered  whenever  similarity  exists. 

In  addition,  the  downstream  boundary  conditions  at  infinity  are 
automatically  satisfied.  The  mapping  also  allows  one  to  obtain  simpler 
boundary  conditions  in  relating  the  optimal  (£,n)  coordinates  to 
conformal  (p,q)  coordinates. 

Since  the  numerical  calculation  is  to  be  carried  out  in  optimal 
coordinates,  the  transformation  relation  can  be  rewritten  in  optimal 
coordinates.  Using  the  coordinate  relation  given  in  Equations  (2.23) 
and  (2.24),  the  transformation  relation  reads  in  optimal  coordinates 
as 
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(4.6) 


and  the  scale  factor  h  is  given  by 

2  — 

h2  *  l[(t}  +  <5t)2]|1+T1 


(4.7) 


It  can  be  shown  that  the  similarity  solutions  are  recovered  in  the  case 
of  a  blunted  wedge  in  terms  of  the  present  variables  £  and  n  if  the 
exponent  n  is  related  to  the  Falkner-Scan  pressure  gradient  parameter 
8  given  by 
2n 
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1+n 


(4.8) 
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It  is  interesting  to  note  that  the  mapping  given  by  the  Equation  (4.6) 
yields  the  limiting  cases  of  the  thin  flat  plate  and  the  vertical 
wall  when  0=0,  Re  =  0  and  8=1  and  Re  ^  0  respectively. 


The  scale  factor  of  the  transformation  in  optimal  coordinates 
is  rewritten  in  nondimens ional  form  as 


h  =  [(c2t2+(^  +  <i w)2)]  2 


(4.9) 


It  is  important  to  note  that  the  independent  variable  n  is  defined 
here  (i.e.  in  Equation  (4.9)),  in  such  a  way  that  n  =  0  corresponds 
to  the  value  of  the  conformal  coordinate  q  at  the  wall  and  it  is 
denoted  by  q^.  With  this  new  definition  of  the  independent  variable  n, 
the  body  is  always  placed  at  n  =  0  in  the  optimal  coordinate  plane 
and  n  =  C  represents  the  displacement  surface  in  optimal  coordinates. 

The  additional  relations  required  for  the  solution  of  the  governing 
differential  Equations  (2.29)  and  (2.30)  are  given  as  follows 
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Along  the  stagnation  streamline,  i.e.  E  »  0,  the  term  h^/Eh  appearing 
in  the  fi^/EH  term  in  the  vorticity  Equation  (2.30)  becomes  undefined 
and  hence,  must  be  evaluated  with  the  proper  limit  process. 


The  limit  process  gives 
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In  order  to  performs  the  boundary- layer  analysis  for  the  blunted 
wedge  oroblem,  the  coordinates  are  normalized  with  respect  to  the 
nose  radius  of  curvature  a,  at  the  stagnation  point,  such  that  the 
wall  location  of  the  blunted  wedge  surface  is  given  as 


(1-6)' 


(4.13) 


This  type  of  normalization  is  adopted  in  order  to  remove  the  explicit 
apperance  of  the  Reynolds  number  Re  from  the  analysis. 

The  scale  factor  of  the  transformation  can  then  be  written  as 
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(4.14) 


and  the  inverse  scale  factor  H  is  given  by 
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Now  the  pressure  gradient  parameter  8^,  appearing  in  the  boundary 
layer  equation  is  given  as 


8 


i 


[1  + 
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(4.16) 


This  completes  the  necessary  details  required  to  solve  the  flow  equations. 

The  cooredinate  equations  required  to  generate  optimal  coordinates 
are  given  by  Equation  (2.25)  and  the  three  boundary  conditions  are 
given  in  Equations  (2.27)  through  (2.28a,b).  The  fourth  boundary 
condition  t(£,0)  is  to  be  determined  iteratively  and  the  iterative 
scheme  required  to  compute  the  missing  boundary  condition  is  given 
in  detail  in  Section  II. 

This  completes  the  mathematical  formulation  of  the  blunted  wedge 
problem  in  optimal  coordinates.  The  details  of  the  numerical  method 
used  to  solve  the  governing  differential  equations  are  presented  in 
Section  III. 


2.  The  Blunted  Plate  Example 

In  this  example  problem,  two  dimensional  symmetric  flow  past  a 
semi- inf inite  blunted  plate  is  obtained  in  optimal  coordinates.  The 
body  geometry  is  shown  in  the  physical  (x,y)  plane  in  Figure  4.  The 
body  is  mapped  to  the  line  q  *  0  in  the  conformal  plane  5 ,  as  shown 
in  Figure  Figure  5,  by  applying  the  Schwarz-Christof fel  transformation 
and  results  in  a  relation  given  by 
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Integration  of  (4.17a)  gives  the  transformation  relation 
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The  x,y  coordinates  are  nondimensionalized  by  the  viscous  length  v/U 
where  v  is  the  kinematic  viscosity  and  is  some  reference  speed. 

Since  the  mapping  is  conformal,  the  scale  factors  in  both  directions 
are  the  same  and  are  given  by 
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which  gives  for  the  present  geometry 
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where 


Re 


L 


2 

IT 


The  Reynolds  number  is  defined  based  on  the  nose  radius  of  curvature 
as  Re  *  U^a/v.  The  nose  radius  of  curvature  at  the  stagnation  point 
can  be  shown  to  be  related  to  the  thickness  of  the  plate  L,  by  the 


following  relation 

<ReL  +  c£3/2 


(4.20) 


where  q^  is  the  coordinate  q  at  the  body  surface.  The  flow  past 
different  body  shapes  can  be  studied  by  a  suitable  combination  of  the 
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parameters  Re^  and  q^  whereas  the  effects  of  Reynolds  number  can  be 
studied  by  varying  Re,  keeping  Re.  and  q  fixed. 

The  conformal  mapping  used  here  transforms  the  flow  in  the 
physical  (x,y)  plane  into  a  stagnation  point  flow  in  the  conformal 
(p,q)  plane.  The  advantages  of  this  type  of  transformation  are  the 
same  as  mentioned  in  the  blunted  wedge  example  problem. 

It  is  interesting  to  note  that  as  a  limiting  case,  as  Re  -*•  0, 

L 

the  transformation  (4.17b)  represents  the  flow  past  a  parabola  which 
corresponds  to  the  flow  past  the  thin  semi-infinite  flat  plate  in  the 
physical  plane.  Therefore  the  analysis  includes  the  flow  past  the 
thin  semi- infinite  flat  plate  as  a  special  case  as  Re^  -»■  0. 

The  governing  differential  equations  and  the  boundary  conditions 
remain  the  same  as  given  in  the  previous  example  problem  by  Equations 
(-.5)  through  (2.81),  except  that  in  order  to  recover  the  correct  . 
asymptotic  form  at  downstream  infinity,  the  value  of  8  in  Equations 
(2.18)  and  (2.44)  is  set  equal  to  zero  to  give  the  Blasius  equation, 
which  is  the  appropriate  downstream  boundary  condition  for  the  blunted 
plate  problem. 

The  scale  factor  h,  in  optimal  coordinates  can  be  written  as 

22  2  1/2 

h(5»n)  -  [(S2t2  +  (qw  +^)  -  ReL)  +  4ReL(qw  +  J)  ]  .  (4.21) 


The  additional  relations  required  for  the  parabolized  vorticity- 
streamfunction  Equations  (2.29)  and  (2.30)  are  given  as 

(4.22) 
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which  can  be  written  in  optimal  coordinates  as 


(4.23) 


Also, 
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As  previously,  we  need  the  following  limit  expression.  Lim  £-*■() 
h^/Ch  can  be  shown  to  be  equal  to 
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For  the  boundary- layer  analysis  the  coordinates  (x,y)  are  normalized 
with  respect  to  the  radius  of  curvature  a.  Then  the  transformation 
relation  (4.17b)  in  nondimens ional  form  becomes 


z  -  \  [C(b  +  S2)1/2  +  b  sin-1  (^j) 
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1/2 

Now  defining  a  =  q^/a  and  substituting  into  the  radius  of 
curvature  expression,  we  can  show  that  the  body  shape  parameter  a,  in 
the  conformal  plane,  is  related  tc  the  thickness  of  the  plate  in  the 
physical  plane  by  the  relation 
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The  shape  of  the  body  for  various  values  of  the  parameter  a,  is  shown 
in  Figure  6.  As  a  special  case  of  the  parameter  a  ■  1,  the  trans¬ 
formation  relation  (4.27)  gives  flow  past  a  parabola. 

The  scale  factor  for  the  boundary-layer  in  optimal  coordinates 
is  written  as 


uV  +  2?V  (2a2  -  a2/3)  +  a4/3]1/2 


(4.30) 


The  governing  differential  equations  for  the  boundary- layer  analysis 
are  the  same  as  given  in  the  previous  example  problem  by  Equations 
(2.37)  through  (2.45)  except  that  in  the  downstream  infinity  condition 
(2.44),  8  is  set  equal  to  zero  to  give  the  appropriate  downstream 
boundary  condition  for  the  blunted  plate  problem. 

The  remaining  relations  required  for  the  boundary-layer  solution 
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As  was  mentioned  in  our  first  example  problem,  the  coordinate 
equations  in  this  case  remain  the  same  as  were  given  by  Equations 
(2.25)  through  (2.28),  except  for  the  surface  boundary  condition  t(£,0). 
Details  of  the  determination  of  the  surface  boundary  condition  t(£,0) 

are  presented  in  Section  II  of  this  analysis. 

This  completes  the  mathematical  modelling  of  the  blunted  plate 
problem.  Details  of  the  solution  procedure  are  discussed  in  the 
section  on  numerical  method  of  solution.  The  results  of  the  present 
analysis  are  discussed  in  the  following  section. 
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V.  RESULTS  AND  DISCUSSION 

In  the  present  study  an  effort  has  been  made  to  demonstrate 
that  a  set  of  optimal  coordinates  dependent  on  the  displacement  effect 
can  be  generated  numerically  for  two-dimensional  symmetric  incom¬ 
pressible  laminar  flow  past  blunt  and  sharp  nosed  bodies.  Various 
critical  aspects  of  the  viscous  flow  analysis  in  optimal  coordinates 
for  high  as  well  as  low  Reynolds  number  were  discussed  in  the  previous 
sections.  The  present  investigation  is  therefore  viewed  as  an  initial 
step  for  generating  flow  dependent  coordinate  systems  numerically. 

The  present  analysis  has  provided  us  with  useful  information  in  order 
to  understand  and  investigate  the  inviscid  and  viscous  interaction 
problem  in  a  different  manner.' 

Present  calculations  were  confined  to  flow  past  semi-infinite 
bodies  under  conditions  such  that  separation  does  not  occur.  The 
reasons  for  choosing  such  body  shapes  are  firstly,  because  they  are 
the  easiest  ones  to  analyze,  flows  involving  separation,  for  example, 
being  much  more  difficult  .  Secondly,  for  comparison  purposes  the 
simple  body  shapes  such  as  a  parabola,  vertical  wall  and  thin  flat 
plate,  possess  previously  obtained  accurate  analytical  or  numerical 
solutions  which  can  be  used  as  a  basis  for  comparison. 

The  solutions  were  obtained  using  the  governing  differential 
equations  developed  for  the  example  problems  in  Section  II  and  numerical 
technique  described  in  Section  III.  The  solution  steps  for  the  coupled 
set  of  flow  equations  and  the  coordinate  equations  are  presented  in 
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Figures  (7-9).  The  results  for  specific  example  cases  are  discussed 
in  the  following. 

1 .  Boundary  Laver  Results  in  Optimal  Coordinates 

la.  Blunted  Wedges: 

In  the  mathematical  formulation  of  the  wedge  problem  for  the 
boundary-layer  analysis  in  Section  IV,  it  was  shown  that  the  flow  past 
different  blunted  wedge  shapes  can  be  characterized  by  a  single 
parameter  8,  related  to  the  slope  of  the  wedge  at  downstream  infinity, 
where  6  is  the  Falkner-Skan  pressure  gradient  parameter.  In  the 
present  analysis,  boundary -layer  solutions  are  obtained  for  3 
between  0  and  1,  where  8  =  1  represents  flow  against  a  vertical  wall, 
commonly  known  as  Hiemenz  flow  and  8*0  represents  the  flow  past  a 
parabolic  cylinder. 

Plots  of  the  pressure  gradient  parameter  8^,  for  various  values 
of  the  downstream  Falkner-Skan  pressure  gradient  parameter  8  are 
presented  in  Figure  10.  Figure  10  shows  that  the  pressure  gradient 
for  the  case  of  the  parabola  approaches  zero  at  downstream  infinity 
and  therefore  the  Blasius  solution  for  flow  past  a  flat  plate  is  the 
correct  solution  asymptotically  at  downstream  infinity.  Figure  10 
also  shows  the  8^  distribution  for  other  values  of  the  downstream 
Falkner-Skan  parameter  8  which  corresponds  to  increasing  the  downstream 
asymptotic  wedge  angle.  If  the  downstream  wedge  angle  is  increased 
to  90®,  the  Hiemenz  flow  or  the  flow  past  a  vertical  wall  is  obtained 
and  the  pressure  gradient  parameter  becomes  constant  over  the  wall 
surface  for  this  flow. 
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The  streamwise  distributions  of  the  skin-friction  coefficient 
gCp**^)  are  shown  in  Figure  11  including  flow  past  the  parabola 
given  by  0  =  0.  For  all  values  of  B,  the  skin  friction  coefficient 
attains  at  the  stagnation  point  a  value  of  1.233  which  is  the  exact 
value  of  skin  friction  coefficient  for  the  stagnation  flow  against 
a  vertical  wall.  For  the  case  of  a  parabola,  the  skin  friction  approaches 
the  flat  plate  solution  at  downstream  infinity.  This  is  as  it  should 
be,  since  the  appropriate  downstream  flow  in  the  case  of  parabola 
is  governed  by  the  Blasius  flat  plate  solution.  With  increasing  B, 
the  skin  friction  coefficient  approaches  the  corresponding  Falkner- 
Skan  solution  at  downstream  infinity.  Observe  that  for  B  *  0.7, 
the  skin  friction  profile  is  very  close  to  that  of  the  flow  against  a 
vertical  wall,  giving  an  indication  that  the  flow  is  approaching  the 
stagnation  flow. 

The  displacement  thickness  functions  for  B  between  0  and  1  are 
obtained  in  optimal  coordinates  and  their  values  in  conformal  coor¬ 
dinates  are  shown  in  Figure  12.  Displacement  thickness  results  obtained 

at  the  stagnation  point  show  very  good  agreement  with  those  of 

15  9 

Davis  as  well  as  those  of  Van  Dyke.  The  displacement  thickness 
for  the  case  of  a  parabola  also  approaches  asymptotically  the  appropriate 
value  obtained  for  the  Blasius  flat  plate  solution  at  downstream 
infinity.  As  expected,  the  displacement  thickness  approaches  a 
constant  value  as  the  downstream  wedge  angle  is  increased  by  increasing 
the  value  of  Falkner-Skan  parameter  B  to  1.  Since  the  displacement 
thickness  in  conformal  coordinates  is  constant,  according  to  optimal 
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coordinate  theory,  the  conformal  coordinates  for  stagnation  flow 
(i.e.  for  6*1),  are  optimal  coordinates. 

lb.  Blunted  Plates: 

It  was  previously  shown  in  Section  IV  (see  Figure  6)  that  for 
boundary-layer  analysis,  the  flow  past  the  blunted  plate  is  characterized 
by  a  parameter  a.  Various  body  shapes  for  the  blunted  plate  problem 
are  shown  in  Figure  6  with  their  corresponding  values  of  the  parameter  a. 
In  the  blunted  plate  case,  a  =  1  represents  a  flow  past  a  parabola 
and  a  =  0  represents  the  flow  past  a  semi-infinite  thick  plate  with 
sharp  corners.  This  type  of  blunt  body  is  considered  in  order  to 
study  the  displacement  flow  past  bodies  with  smooth  as  well  as  blunted 
shoulders  and  to  investigate  the  effects  of  such  convex  corner  flows 
on  the  optimal  coordinate  configuration. 

The  surface  skin  friction  distributions  shown  in  Figure  13  are 
presented  for  the  bluntness  parameter  a  between  0.2  and  1.  The  skin 
friction  variation  for  flow  past  the  parabola  has  already  been  discussed 
in  the  case  of  the  blunted  wedge  problem  and  the  results  of  the  same 
are  presented  here  for  the  sake  of  comparison  with  the  skin  friction 
variation  for  different  values  of  the  parameter  a.  Figure  6  shows 
that  decreasing  the  value  of  the  parameter  a,  results  in  a  body  shape 
having  larger  curvature  around  the  corner.  The  viscous  flow  past  such 
bodies  would  experience  larger  and  larger  turning  angles  around  the 
corner.  This  fact  is  displayed  in  the  skin  friction  plots,  which  show 
that  with  decreasing  values  of  a,  the  skin  friction  shows  an  earlier 
drop,  larger  gradients,  and  a  smaller  region  of  adjustment  to  the 


downstream  infinity  value  as  compared  to  the  skin  friction  variation 

for  flow  past  a  parabolic  surface.  Similar  features  are  also  displayed 

by  other  flow  quantities  such  as  pressure  gradient  and  displacement 

thickness  for  decreasing  values  of  o,  in  Figures  14  and  15  respectively. 

These  characterstic  features  of  the  flow  quantities  such  as  are  shown 

in  the  skin  friction  and  the  displacement  thickness  are  consistent 

with  che  pressure  gradient  plots  obtained  for  the  viscous  flow  past 

the  corresponding  body  shapes,  represented  by  the  parameter  a.  One 

of  the  flow  cases  which  we  will  study  further  is  the  flow  past  a 

blunt  body  represented  by  a  =  0.2  which  will  turn  out  to  be  a  case 

near  separation.  A  classical  boundary-layer  analysis  of  flow  past  a 

body  which  is  represented  by  a  less  than  approximately  0.2  would 

indicate  the  occurance  of  flow  separation  aft  of  the  corner.  This 
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fact  has  previously  been  observed  by  Davis  in  his  interacting  flow 
analysis  for  high  Reynolds  number  flow  past  a  body  represented  by 
a  ■  0.1.  Considerable  numerical  difficulties  were  encountered  in  order 
to  obtain  converged  solution  for  the  flow  past  a  blunt  body  given  by 
a  less  than  0.2.  The  pressure  gradient  parameter  in  such  cases  reaches 
limiting  values  below  -0.199  and  the  shear  stress  drops  to  zero  with 
an  Infinite  slope  aft  of  the  corner  indicating  the  onset  of  separation 
at  that  point  and  no  solution  could  be  achieved  beyond  the  separation 
point  indicating  that  the  present  solution  technique  must  be  modified 
if  one  wishes  to  compute  separated  flows  correctly.  We  will  see  in  the 
later  part  of  our  discussion  that  the  optimal  coordinates  for  the  case 
of  a  ■  0.2  exhibit  critical  behavior  for  this  near  separated  flow  case. 


lc.  Optimal  Coordinates  for  Blunted  Wedges  and  Blunted  Plates: 

It  was  mentioned  in  the  preceeding  section  that  the  displacement 
thickness  is  constant  for  flow  past  a  vertical  wall  (Hiemenz  flow). 

The  condition  for  optimal  coordinates  requires  that  the  displacement 
thickness  should  be  constant  in  optimal  coordinates.  This  implies 
that  the  conformal  coordinate  system  should  be  optimal  for  the 
stagnation  flow.  The  validity  of  the  statement  is  justified  from 
the  fact  that  the  coordinate  equation  for  optimal  coordinates  yields 
a  constant  solution  equal  to  1  for  the  coordinate  function  t(£,n) 
in  the  entire  flow  field. 

The  optimal  coordinate  (£,n)  system  and  the  conformal  (p,q) 
coordinate  system  for  values  of  the  wedge  parameter  8  between  0.7 
and  0.0  are  presented  in  Figures  16  through  18.  The  optimal  coordinates 
for  the  wedge  parameter  8  =  0.7  show  little  difference  from  the 
conformal  coordinate  system  as  shown  in  Figure  16a  because  the  flow 
represented  by  8  *  0.7  has  characteristics  very  close  to  stagnation 
flow.  The  differences  between  the  two  coordinate  systems  disappear 
as  we  go  far  from  the  body  which  is  obvious  from  the  coordinate  plots 
of  Figure  16b  for  flow  past  the  same  body  shape.  With  a  decrease  of 
the  wedge  angle  at  downstream  infinity,  thereby  decreasing  the  para¬ 
meter  8,  the  initial  differences  between  the  optimal  coordinates  and 
the  conformal  coordinates  near  the  surface  of  the  body  increase  and 
this  difference  dies  out  as  we  go  far  from  the  body.  This  can  be 
seen  in  the  optimal  coordinate  plots  shown  in  Figures  17a  and  17b  for 
the  flow  given  by  the  parameter  8  ■  0.5.  Another  interesting  feature 


63 


displayed  in  the  optimal  coordinate  plots  for  8  =  0.5  in  Figure  17a 
is  that  the  variations  in  the  n  direction  are  larger  than  the 
variations  in  the  £  direction.  Similar  features  are  also  exhibited 
in  the  optimal  coordinate  plots  for  8  =  0.3  as  shown  in  Figures  18a 
and  18b  except  that  the  optimal  coordinates  show  larger  initial 
differences  with  respect  to  the  conformal  coordinates  near  the  body 
surface.  This  difference  subsequently  disappears  as  the  downstream 
conditions  are  reached.  The  optimal  coordinates  for  the  flow  past 
the  parabola  are  presented  in  Figures  19a  and  19b.  The  optimal 
coordinates  approach  the  corresponding  conformal  coordinate  system 
asymptotically  at  downstream  infinity.  This  characteristic  of  the 
coordinate  function  is  in  complete  agreement  with  the  asymptotic 
approach  of  the  downstream  condition  of  the  displacement  flow  past  a 
parabola. 

Next,  the  optimal  coordinates  are  obtained  for  blunted  plates 
for  various  values  of  the  parameter  a.  Typical  optimal  coordinates 
are  presented  for  a  in  the  range  between  0.2  and  0.7  in  Figures  20 
through  23.  The  plots  of  the  coordinates  near  the  surface  of  the 
body  are  shown  in  Figures  20a,  21a,  22a,  and  23a.  There  are  larger 
differences  between  the  optimal  coordinates  and  the  corresponding 
conformal  coordinates  near  the  surface.  The  differences  between  the 
two  coordinates,  as  anticipated,  disappear  near  downstream  infinity. 
This  is  clear  from  the  plots  of  the  coordinates  shown  in  Figures  20b, 
21b,  22b,  and  23b.  Also  the  gradients  of  the  curve  represented  by 
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5  *  C  (where  C  is  a  constant) ,  near  the  surface  and  around  the 

shoulder  region  decrease  with  an  increase  of  shoulder  bluntness  of  i 

the  blunt  body. 

2a.  Parabolized  Navier-Stokes  Results  in  Optimal  Coordinates: 

The  parabolized  vorticity-stream  function  results  are  obtained 

/ 

for  the  blunted  plate  and  the  blunted  wedge  problems  for  low  as  well  as  high 
Reynolds  number  cases.  The  numerical  solutions  are  obtained  with  a 

step  sizes  of  AS  =  1/40  and  AN  =  1/100.  For  a  few  cases  of  the  * 

problem  parameters,  the  solutions  are  obtained  with  larger  step 

sizes  of  AS  =*  1/20  and  AN  *  1/60.  In  order  to  solve  the  coupled  set 

of  coordinate  and  flow  equations  two  different  approaches  of  treating 

the  partial  differential  equations  were  explored.  In  one  case 

the  solution  of  the  parabolic  vorticity-stream  function  equations, 

and  the  *  (star)  level  step  of  the  ADI  scheme  to  solve  the  coordinate 

equations  are  carried  out  simultaneously.  At  each  station  represented 

by  5  *  constant,  the  condition  of  the  optimal  coordinates  is  enforced 

before  advancing  to  the  next  station  in  the  streamwise  direction. 

During  this  operation  the  (n+1)  level  coordinate  distribution  is  kept 
unaltered.  Once  the  condition  of  optimal  coordinates  is  satisfied  by 
the  edge  quantities  of  the  flow  equation  in  the  entire  flow  field, 
the  (n+1)  level  step  of  the  ADI  scheme  is  performed.  The  above  steps 
are  repeated  to  obtain  the  converged  solution  of  the  coordinate 
equation. 

In  the  second  approach  the  flow  equations  are  uncoupled  from  the 
coordinate  equation  Iteration  by  linearizing  the  coefficients  dependent 
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on  the  coordinate  function  t(£,n)  by  using  the  preceding  iteration 
value  of  the  coordinate  function.  The  condition  of  optimal  coordinates 
is  not  checked  until  the  coordinate  solutions  are  iterated  to  give 
a  converged  solution.  The  flow  charts  of  both  the  numerical  approaches 
are  presented  in  Figures  8  and  9. 

The  first  scheme  described  here  was  very  successful  in  obtaining 
accurate  converged  solution  for  the  blunted  wedge  cases.  Considerable 
difficulty  was  encountered  with  this  scheme  in  obtaining  solutions 
for  the  blunted  plate  problem.  However,  the  second  scheme  was 
successful  in  achieving  converged  solutions  for  the  blunted  plate  cases. 
In  general,  both  schemes  displayed  poor  convergence  on  the  coordinate 
equations  for  high  Reynolds  number  flow  cases,  especially  in  the  case 
of  the  blunted  plate  problem. 

The  fictitious  time  step  Av  used  in  the  ADI  scheme  was  found  to 
vary  with  the  body  shape  parameter  8  or  a  and  the  Reynolds  number  Re. 

A  time  step  of  10  was  found  to  be  optimum  for  obtaining  a  converged 
solution  of  the  problems  for  low  and  moderately  high  Reynolds  numbers 
of  up  to  100.  Beyond  this  value  of  Reynolds  number  a  larger  time 

4 

step  was  necessary.  For  example,  for  Reynolds  number  equal  to  10  ,  a 
time  step  of  100  was  used  to  obtain  the  converged  solutions.  It  was 
also  observed  that  the  convergence  process  was  very  much  improved  by 
taking  initially  a  larger  time  step  and  subsequently  reducing  it  as 
the  iteration  proceeded.  A  further  improvement  in  convergence  was 
achieved  by  suitable  averaging  of  the  two  levels  of  solutions  and  using 
the  averaged  values  of  the  solution  vectors  at  the  beginning  of  the  next 


iteration. 


2b.  Blunted  Wedges 


The  streamwise  distributions  of  the  vorticity  function  gfp.q^)  at 

the  wall  for  the  case  of  the  parabola  are  presented  in  Figure  24  for 

4 

Reynolds  numbers  between  0  and  10  .  As  the  Reynolds  number  increases, 
the  value  of  the  vorticity  function  gCp.q^  at  the  stagnation  point 
approaches  the  stagnation  flow  value  given  by  the  Hiemenz  solution.  At 
Re  *  10^  the  skin  friction  value  approaches  the  stagnation  flow  value 
of  1.233  and  the  Navier-Stokes  solution  matches  with  the  present 
boundary- layer  solution  at  the  stagnation  point  for  the  flow  past  a 
parabola.  The  variations  of  the  vorticity  function  for  the  thin  flat 
plate  case,  which  is  obtained  from  the  flow  past  a  parabola  with 
Re  »  0  is  also  shown  in  Figure  24.  A  comparison  between  the  present 
results  and  those  of  Davis^  for  leading  edge  skin  friction  on  a 

4 

parabola  is  given  in  Figure  25  for  Reynolds  numbers  between  0.1  and  10  . 
The  agreement  for  all  Reynolds  numbers  is  good.  There  are  small 
differences  between  the  two  results  at  moderately  high  Reynolds  number 
between  300  and  1000  as  shown  in  Table  1.  Figures  26  and  27  show  the 
vorticity  function  distribution  along  the  surface  of  the  blunted 
wedges  for  B  *  0.7  and  0.5  respectively.  Results  of  the  present 
calculations  at  the  stagnation  point  for  B  ■  0.7  and  0.5  agree  well 
with  the  corresponding  results  of  Davis,  U.  Ghia  and  K.  Ghia 
obtained  by  solving  the  full  Navier-Stokes  equations.  At  downstream 
infinity  the  results  of  the  present  calculations  approach  the  boundary- 
layer  limit  exactly. 
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2c.  Blunted  Plates : 


The  blunted  plate  analysis  using  parabolized  Navier-Stokes  equa¬ 
tions  is  characterized  by  two  parameters,  the  bluntness  parameter  a 
and  the  Reynolds  number  Re  based  on  the  nose  radius  of  curvature  at 
the  stagnation  point.  Figure  28  shows  typical  results  for  skin 
friction  for  the  blunt  plate  given  by  a  =  0.7.  Retaining  the  value 
of  a  and  increasing  the  Reynolds  number  Re,  the  shoulder  influence  is 
further  removed  from  the  point  (0,q^).  The  leading  edge  will  appear 
more  like  a  vertical  surface  for  the  oncoming  stream  for  high  Reynolds 
number  flow  cases  and  the  flow  in  the  vicinity  of  the  point  (O.q^) 
will  approach  the  stagnation  flow.  Therefore,  it  is  found  that  with 
an  increase  of  the  value  of  Re,  the  vorticity  value  at  the  point  (O.q^) 
approaches  the  stagnation  flow  value.  The  vorticity  function  g(p,qw) 
drops  along  the  surface  from  the  stagnation  point  value  to  the 
Blasius  value  for  a  flat  plate  at  downstream  infinity.  Decreasing 
the  value  of  the  parameter  a  decreases  the  bluntness  of  the  shoulder. 

The  skin  friction  distribution  for  a  =  0.5  as  shown  in  Figure  29 
displays  larger  gradients  and  an  earlier  drop  from  the  stagnation  point 
value  as  compared  to  the  skin  friction  distribution  for  the  parameter 
a  *  0.7  for  the  same  Reynolds  number. 

2d.  Optimal  Coordinates  for  Blunted  Wedges  and  Blunted  Plates: 

In  Che  case  of  boundary- layer  analysis  we  have  analyzed  the  effects 
‘  'he  body  shape  parameter  on  the  optimal  coordinates.  In  this  case 
►  -"•<  ,*  of  varying  the  Reynolds  number  Re  have  been  investigated. 

'  ■  «■  o<<rd* nates  for  flow  past  blunted  wedges  and  blunted  plates 
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have  been  obtained  for  laminar  unseparated  flow  cases.  The  coordinate 
calculations  were  specialized  to  the  limiting  cases  of  flow  past 
a  parabolic  cylinder,  the  thin  flat  plate  and  the  flow  against  the 
vertical  wall. 

The  streamwise  distributions  of  the  coordinate  function  for  the 

case  of  a  parabola  are  shown  in  Figure  30a  for  various  Reynolds 

4 

numbers  Re,  which  vary  from  0  to  10  .  The  coordinate  function  t(C,0) 

monotonically  approaches  the  correct  downstream  condition  at  infinity. 

With  increasing  Re,  the  variations  approach  the  boundary-layer 

distribution  and  results  of  t(£,0)  match  exactly  with  the  boundary- 

4 

layer  solution  at  Re  =  10  as  shown  in  Figure  30b,  except  near  the 
leading  edge.  For  the  case  of  a  thin  flat  plate,  the  optimal 
coordinates  are  obtained  by  setting  Re  =  0  in  the  parabola  problem. 
Results  of  the  surface  distribution  of  t(?,n)  for  the  thin  flat  plate 
are  also  shown  along  with  the  results  for  the  parabola  in  Figure  30a. 
Effects  of  Re  on  the  optimal  coordinates  for  increasing  value  of  the 
downstream  wedge  angle  are  presented  in  Figures  31  and  32  for  wedge 
parameters  of  8  ■  0.3  and  0.5,  respectively.  Similar  features  as 
described  in  the  case  of  the  parabola  are  observed  in  these  cases. 

The  coordinate  distribution  for  6  =  0.7,  which  represents  a  flow 
closer  to  a  stagnation  flow,  is  shown  in  Figure  33.  As  anticipated, 
the  parabolized  Navier-Stokes  solutions  agree  with  the  boundary-layer 
solutions  at  lower  Reynolds  number  than  the  cases  discussed  earlier. 

The  optimal  coordinates  were  obtained  for  the  blunted  plate  cases 
for  Reynolds  number  between  0.1  and  100.  Reynolds  numbers  larger  than 


a  value  of  100  were  not  obtained  because  of  poor  convergence  of  the 
present  numerical  scheme.  Typical  plots  of  the  surface  value  of  the 
coordinate  function  are  shown  in  Figures  34  and  35  for  a  =  0.7  and 
0.5  and  for  Reynolds  numbers  between  0.1  and  100.  In  all  these  cases 
with  increasing  values  of  Re,  the  coordinate  function  at  the  surface 
tends  to  assume  a  boundary-layer  variation.  Since  a  =  1  represents 
the  case  of  a  parabola,  the  optimal  coordinate  results  for  this  case 
have  already  been  discussed  along  with  the  blunted  wedge  cases  and 
will  not  be  repeated  here. 

The  conformal  and  the  optimal  coordinate  plots  for  the  thin  flat  plate 
and  the  parabola  are  shown  in  Figures  36(a,b)  and  37 (a, b)  respectively. 
The  initial  differences  between  the  two  coordinate  systems,  such  as 
are  observed  in  Figure  36a  for  the  flat  plate  case  and  in  Figure  37a 
for  the  case  of  a  parabola,  tend  to  reduce  as  we  go  further  downstream 
from  the  body  surface.  This  is  clear  from  the  plots  shown  in  Figures- 
36b  and  37b. 

Figure  38a  shows  the  conformal  and  the  optimal  coordinates  near  the 
surface  of  a  blunted  wedge  represented  by  6  ■  0.5  and  for  Reynolds 
number  Re  ■  100.  Further  downstream  of  the  blunted  wedge  surface, 
the  optimal  coordinates, . as  shown  in  Figure  3Sb,  tend  to  merge  with 
the  corresponding  conformal  coordinates. 


VI.  CONCLUSIONS 


In  the  present  analysis  of  viscous  flow,  it  has  been  successfully 
demonstrated  that  it  is  possible  to  numerically  generate  a  particular 
set  of  flow  dependent  coordinates  called  optimal  coordinates.  Results 
of  the  integration  of  these  coordinate  equations  reveal  that  the 
optimal  coordinates  used  here  may  be  viewed  as  a  generalized  form  of 
conformal  coordinates.  The  optimal  coordinates  considered  here 
automatically  degenerate  into  convenient  sets  of  coordinates  for 
specific  flow  fields;  for  example,  in  the  case  of  first-order  boundary- 
layer  flow  past  a  thin  flat  plate,  the  optimal  coordinates  degenerate 
into  parabolic  coordinates,  which  are  the  most  natural  form  of  the 
coordinates  for  the  thin  flat  plate.  Similarly  for  the  case  of  flow 
against  a  vertical  wall,  the  optimal  coordinates  become  simple 
Cartesian  coordinates,  for  both,  the  boundary-layer  flow  and  the 
Navier-Stokes  flow. 

The  advantages  of  viscous  flow  analysis  in  optimal  coordinates 
are  as  follows.  Firstly,  the  boundary-layer  solution  obtained  in 
optimal  coordinates  is  uniformly  valid  throughout  the  entire  flowfield. 
Secondly,  similarity  conditions  are  automatically  satisfied  when 
similarity  exists,  with  the  present  form  of  the  coordinates. 

The  results  of  the  optimal  coordinate  calculations  using  parabolized 
Navier-Stokes  equations  for  the  case  of  the  parabola  and  blunted  wedges 
have  shown  excellent  agreement  with  the  existing  results  using  full 
Navier-Stokes  equations  for  all  Reynolds  numbers.  This  indicates 


that  the  parabolic  approximation  to  the  Navier-Stokes  equations 

is  a  more  accurate  representation  of  the  full  Navier-Stokes  equations 

if  the  equations  are  generated  in  optimal  coordinates.  The  works 
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of  Davis  and  U.  Ghia  and  Davis  which  used  the  parabolized  Navier- 
Stokes  equations  in  conformal  coordinates  (which  are  semi-optimal 
in  nature)  have  shown  reasonable  agreement  with  full  Navier-Stokes 
results  even  for  low  Reynolds  numbers,  and  therefore  it  was  to  be 
expected  that  the  use  of  optimal  coordinates  in  the  parabolic  form 
of  the  Navier-Stokes  equations  would  further  improve  the  results. 

The  results  of  the  present  calculations  verify  this  assumption. 

The  governing  differential  equations  and  the  numerical  scheme 
developed  in  the  present  study  for  boundary- layer  analysis  was  found 
to  be  successful  in  predicting  first  order  optimal  coordinates. 

Present  numerical  calculations  showed  that  the  cross-derivative 
term  ,  in  the  coordinate  equation  can  be  treated  explicity  or 
implicitly  in  the  numerical  solution  of  the  coordinate  equations  for 
boundary-layer  analysis,  whereas  it  was  found  to  be  essential  to 
treat  the  same  cross-derivative  term  8^  implicitly  for  the  parabolized 
Navier-Stokes  analysis  for  the  low  Reynolds  number  cases. 

In  the  parabolized  Navier-Stokes  analysis,  the  coordinate 
equation  solution  scheme  encountered  certain  difficulties  which  are 
worth  mentioning.  The  first  iterative  scheme  mentioned  in  the  parabolic 
Navier-Stokes  analysis  displayed  poor  convergence  of  the  solution  of 
the  coordinate  equation  for  high  Reynolds  numbers.  However  the  use  of 
the  second  or  the  uncoupled  iterative  scheme  was  found  to  improve 


the  convergence  of  the  problem  for  most  blunt  body  cases.  It  was 
also  observed  in  general  that  the  relaxation  process  of  the  coordinate 
equation  becomes  very  sensitive  to  the  fictitious  time  step,  Av, 
and  the  initialization  process  of  the  coordinate  function  for  high 
Reynolds  number  cases.  While  the  second  scheme  enjoys  a  greater 
flexibility  of  choice  of  the  time  step  Av  and  the  initialization 
process  of  the  coordinate  function  as  compared  to  the  first,  diffi¬ 
culties  were  still  encountered  for  Reynolds  numbers  greater  than 
100  for  blunted  plate  problems.  Therefore  a  future  effort  should 
be  made  to  resolve  the  problem  of  poor  convergence  for  high  Reynolds 
number. 
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PIG.  8.  FLOW  CHART  FOR  SOLVING  PARABOLIZED  NAVIER-STOKES 

EQUATIONS  IN  OPTIMAL  COORDINATES  BY  COUPLED  SCHEME 
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FIG.  16a.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 


FIG.  16b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 


17b.  CONFORMAL  AMD  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 
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FIG.  18a.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 


FIG.  18b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 
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FIG.  19a.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  PARABOLA 
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FIG.  19b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  PARABOLA 
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Y-AXIS 

FIG.  20b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


PIG.  21b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


FIG.  22a.  CONFORMAI.  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


FIG.  22b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 
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FIG.  23a.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


CONFORMAI.  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  PLATE 


31.  OPTIMAL  COORDINATE  FUNCTION  DISTRIBUTIONS  ON  BLUNTED  WEDGE  SURFACE 
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FIG.  37b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  A  PARABOLA 


FIG.  38a.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR  BLUNTED  WEDGE 
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FIG.  38b.  CONFORMAL  AND  OPTIMAL  COORDINATES  FOR 
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